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TROL ASPECTS OF CAUCA VALLEY DEVELOPMENT 
Carlos S. Ospina,“ F. ASCE 


mes, 
The Cauca Valley project applies two techniques relatively new in Latin 
America: integrated water - resource development and administration by a 
decentralized, TVA-type agency. The flood-control aspects are greatly af- 
fected by the other aspects, mainly power, irrigation anddrainage. 
@ Parallel matters discussed herein are the hydrology of the region, which | . 
1 is in the equatorial “doldrum belt, ¥ ’ and the rate. of ——— as affected = 
a. by economic ic and political considerations. is 


INTRODUCTION 


lanked by the two acetic ranges of the » Colombian 1 Andes, the Cauca 
Valley is a flat, verdant strip of land 200 km long, an average of | 20 km wide 
and almost 400,000 ha in extent; see Figs. 1 and 2. " (Abbreviations and equ equi- 
valents are listed atthe endofthis paper.) = 
: As the elevation is 950 m above sea level, it enjoys a pleasant, balmy a - 
a climate, despite the equatorial latitude. Long- -term averages are: 
24°C, humidity 63%, and annual rainfall 1,200 mm. 
: Although well settled for over four centuries, only a third of the valley is = 


cultivated at present. This low degree of development results fromacom- | 4 
_ plex of sociologic and eer reasons, some of which are discussed in es 4 
‘The Corporacion Auténoma Regional del Cauca (Autonomous 


: Discussion open until ieiaes. 1, 1960. To extend the closing date one month, 
i written request must be filed with the Executive Secretary, ASCE. Paper 2147is 
part of the copyrighted Journal of the Hydraulics Division, Proceedings of the 4 
American Society of Engineers, Vol. 85, No. HY 9, September, (1959, Meo 


1, Associate, Tippetts- -Abbett- Stratton and 
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FIGURE 1. VICINITY MAP 

David E. Lilienthal, former chairman of the TVA, assisted the government in es 
preparing the basic legislation; for this reason, the Develop- 


resource developments are planned within three 

q (provinces): Cauca, Valle and Caldas: Hence, the abbreviation “CVC” 
_ common local usage. The CVC is ound with the integrated =<. 

a of water resources including flood control, electric power, drainage, , irriga- i. 

tion, pollution abatement and related uses. The CVC in 1955 engaged the a 

consulting engineering firms of Olarte, Ospina, Arias & Payan, Ltda. of _ a, 


Bogota, Colombia, and Gibbs & Hill, Inc., and Tippetts- - Abbett - McCarthy, both 
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“quent flooding has discouraged significant developments within | the extensive 
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of New York, to carry out the planning and engineering required for the Cauca . 


: _ Valley development. ~ This paper is based on the studies made by these three ie 
together with the technical staffofthe CVC. 


‘ere 


- Within CVC’s development area reside 2,600,000 persons. The center © of . 
the valley is Cali, a booming city of 400, 000, with 2 an annual population in- 
crease of 8%. A hundred kilometers to the west is the Pacific ect al Port 
: Almost a third of the valley is in the flood zone. river 
- most every year, in some years as many as three or four times. A major — E 
_ flood occurs once in 7 to 10 years when about 85,000 ha are inundated by th the 4 
ri (see Figs. and 4), plus about 30,000 ha by tributary streams. 
_ Because the banks of the Cauca River are higher than the adjoining land on 
‘each side, , damages from Cauca River floods extend well beyond the 85,000- ha 
area subject to direct flooding in two ways: me) Floodwater trapped behind the 
river banks causes high groundwater levels over extensive adjacent areas (a 


outlet conditions, flooding extends over a ™ area at the mouths of many of — 


Despite the area flooded, actual is small because the fre- 


outlets, even in the absence of general flood control). — (2): Owing to unfavorable 


_ low-lying areas along the river. The early Spanish settlers placed » 
towns on high ground d for sanitary reasons. . The present u use of land in the 
-_ ‘flood zone and adjacent 2 areas is limited to cattle grazing; when floods occur | 
the cattle are simply moved to higher ground and brought back after the land # 
ee out. Thus most of the benefits from flood protection are those that will 
- result in the future from intensification in use of the land (for agriculture, in- 
dustrial sites and for urbanization) rather than from avoidance of damage to — 
‘Natural valley storage is high because of the large area flooded. During -¥ 
_ the February 1950 flood, 31,000 ha upstream of Cali were submerged beneath = 
‘million ms of water, the average depth being 2. 4 m. n. The siderably 


It is evident that ‘flood control in upper er portion above Cali), 


Characteristics 


Meteorology 


Table 1 gives average at three representative station 
_ Popayan and Piendamé are in the “upper basin ” (that portion of the drainage Ei 
area upstream of the Timba | damsite at La | Balsa), while La Manuelita — 


representative of the valley. ‘The valley | has two distinct dry seasons; the up- 
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16.9 313 14. 8 11.4, 


853 1 100.0 2,121 100.0 1,120 100.0 


a _ Storm rainfall in the upper Cauca basin is spotty and erratic. There are “a 


| _ variations with respect to time and place even during short periods of 5 


less than a day. The precipitation records show that hat heavy rain: rainfalls at the © a 


Another characteristic of storm rainfall in the upper Cauca basin a its n¥ 


| :. low depth as compared With the annual rainfall. Thus, for the 18 years of — / 


6% of the average annual rainfall. This low percentage can be contrasted with 
the much higher percentages observed in other parts of the world, outside the - 
_ latitudes of the equatorial trough (or “doldrums”)in which the upper Cauca 3 
basin is located. In these latitudes, winds are mild and tropical hurricanes 
= are unknown. (Tropical hurricanes universally originate at about 10° to Aa a, 
south latitude and move ‘respectively | northward and southward. ) Moreover, 
3 cold fronts, s such as those that cause ‘major storms in the Canal Zone, do not - 
: reach these latitudes as they never penetrate further south than about 8°N ~- 


ee records at Popayan, the maximum 24- hr rainfall has never exceeded 


‘The available records for the upper basin show that depths of storm rain- 
_ fall are surprisingly low and do not increase with frequency as much as might ea 
_ be expected; see Table 2. The 10-year 24-hr rainfall at Popayan of 100 mm ee 


— about as D. C the 100- -year 24 at 
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Popayan (estimated from an 18-year record) is niall less at 

Bee oy the corresponding basin rainfalls are only about 40% of the point 

rainfall at Popayan. _ Admittedly, the Popayan record is not long enough for a 5 
_ firm conclusion; | on the. ‘other hand, t that the infrequent, » Say 100-year, storms» 
; are not much worse than the 10- -year storms is confirmed from the records —— 

_ Of floods, described below, which indicate that since 1916 all major floods _— 
a _ have be been of about the s same order of ‘magnitude; nor is there evidence of any = 
ear 
Valley inthe sixteenth century, = 

A long-term record (1900 to date) for La Manuelita indicates a rough cor- “a F 
~ relation between annual rainfall in the valley as against the upper basin; this — - 
record also indicates that there is a 6 to 7-year cycle of unusual wet and un- — a 
apparent explanation for this is a periodic shifting of the equa-_ 
trough of low-pressure (inter-tropical convergence or “ITC”)andof 
the accompanying pattern of atmospheric circulation off the Pacific coast of ae 
- the northern part of South America. There is a parallel es of — 
general area where the cold northern-flowing Humboldt current meets a acl 

| = warm current coming from Panama; both | currents then move \ westward to al 
the Galapagos Islands. See Fig. 1. 

bc On the basis of gage records (available on. 

_ in 10 years. - Preceding the beginning of gage records, there were four a 

“historic floods”: in 1916, 1932, 1934 and 1938. The available evidence indi- 

Bs cates that at least three of these were about equal (in peak discharge) t to the 


flood of 1950, confirming a flood « of of the 
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February 1950 has a frequency 10 least ‘with respect to 
‘The worst flood of record at La Balsa had a discharge of 1,225 m3/s from 
a drainage area of 5 »480 km?, The ‘corresponding Myers ‘rating C, in the _ 
_ Jarvis-Myers formula Q=CVA » with Q and A in English units (cubic feet i “ 
per second and square miles), is is 1 ,050. The available records indicate that — 
this value applies to the tributary streams as well. 


values observed in the United States, C is low in the Cauca 


is practically non- ~existent in the Cauca Grainage basin, 
by: a b combination of rain and snow- melt are unknown. 


a When the sane area consists of cities or towns, a high degree of p pro- 2A 
tection is justified to avoid catastrophic leases that may result from great 
floods once in say 100 to 200 years. Ih the of the Cauca Val 
F low, only agricultural land is involved. — Following protection, the land in the — 
flood zone will be devoted to crops (such as rice, corn, beans and cotton) and 7 
to intensive pasture, for which the average growing season (including: harvest- 


_ Sugarcane, for which the growing season is 15 to 18 months, is widely bai 
grown in the Cauca Valley; harvesting takes place the year round—an impor- & 
tant advantage in that the refineries can operate almost continuously —and its = 
_ cultivation will surely be extended. However, it is estimated that not more | ae 
- than 10 to 20 percent of the flood zone will be devoted to cane in the future: Ripa 
_ first, because the heavier soils located therein are better suited to other bo dfeell 
crops, mainly rice for which the high irrigation water- -requirement will be 1s 
relatively easy to supply in the flood zone; second, ample area, mostly free ios 
_ of flood risk, is available for cane cultivation outside the flood zone. : 
Thus, for the expected cropping pattern, 6 months may be assumed to be pe 
_ “te average growing season. . If protection is provided against a flood with 
c> average frequency of once in 10 years (once in 20 growing seasons), 19 crops 
“2 out of 20 will be protected or 95% of the crops; if against a 20-year flood, 39 
2 out of 40 or 97.5% of the crops. ‘Thus, for a flood of half the frequency with ; 


— corresponding expense of protection on the ret of (50% greater), the no 

The foregoing assumes that the area and the crops damaged 


by the 20- ~year flood, would be about the same as by the 10-year flood. ee 


is generally true for the largest part of the flooded area, located upstream and a 
" Following the reasoning | presented above, protection against a 5- -year flood — 


‘would constitute 90% protection so that the benefit gain of 10-year over 5-year | 
4 _ protection is s theoretically only 5%; however, itis necessary t to consider what 3 
- effect the degree of risk will have on the attitude of farmers. Conversations ~ 
; with them have indicated that 10-year protection corresponding to 5% risk leet 
| a (one in twenty) is sufficiently great so that they will not be discouraged from 
Planting crops, but that 5- -year prctection providing rie} risk (one in ten) ody 
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that that obtained by providing control of the 10-year flood. 
There are no important towns or villages within the Cauca River wr flood zone. 
_ Some towns do suffer floods from tributary streams but these involve local | 
projects generally requiring only minor works to achieve protection. ; 
city of Cali, although not suffering serious flood damage at present, 
: expanded into the flood plain to the extent that further expansion has been : 
limited to narrow flood-free zones to the north and south; see Fig. 9. The 
_—Aguablanea project, described below, is being constructed to protect 5,000 ha aj 
i of land adjacent to the city against a flood having an estimated frequency of ; 
100 years, under present conditions. Following completion of the Timba _ <4 
% reservoir and the Cauca River channel improvement, the degree of protection ire 
will be against a flood of at least 200-year 


"Methods of Obtaining Flood. Control = 
a o a -eam of the upp upper —_ of the valley at La a Balsa, t there are two reser- 
voir sites: : Timba and Salvajina. The drainage areas controlled are, respec 


_ 


> 


tively, 5,480 and 3,960 km2. Other significant areas sare: 
At Juanchito near Cali 9, 060 km* 


The ‘Timba reservoir is more valuable as a flood-control structure be- 
cause it controls 38% more drainage area than Salvajina. Moreover, various 
_ additional factors » which cannot be discussed here because of space limitations, — 

- have required that, of the two ‘reservoirs, Timba be given a higher priority. nit 

_ The control of floods on the tributaries must be erat 3 jointly vonage the i 


a as discussed i in a later section. — 


i ‘The desired degree « of flood protection is obtainable by storage in the Ti Timba 
‘reservoir and improvement of the channel of the Cauca River, the optimum = 
a combination being that giving the lowest overall cost, considering the yates of of 
i co arn nner energy | lost because of storage reserved for flood | control. _— 
i a — channel cap capacity of the river is about 550 m3/s at La Balsa and — 
650 m3/ s at Cali. The amount of reservoir storage necessary to reduce past — 
floods | of record down to Prnepataion limits was determined, and this study resulted 
in the e development of a ‘rule curve,” ‘similar to that shown in Fig. 5, giving 
the reservoir level which may not be exceeded on a given calendar date. sei 


two major hydroelectric stations (Anchicaya and Calima—See Fig. 2) and — 
- coal -burning t thermal station with ani installed i capacity of 53 000 kw kw located & 

at Yumbo, a suburb of Cali. From a study of system power-output, using a as) om 
basic data streamflow records from 1945 to date, it was determined that aa >. 


firm power "000 kw. ‘The t from Timba, if the rule curve in Fig. 5 is oe 
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4 
er ruie Curves Were develop or lesser degrees Of reductio 
_ Besides flood control and power, Timba will provide the means Oi divert- 
ing water to two principal irrigation canals, one on each side 


‘Reservoir level Action required 
Higher: than Rule Curve Moke the total outflow equal 
lower thon E1.1025.0 the _ lowest of the following two 
Balso and Coli 
Open all outlets to prevent ao 
ferther rise in reservoir level 
E~ 
-FIGU RE 5. 5. RULE CURVE, TIM BAR 
‘The: total net area in the valley to 280,000 ha, of which 
q 113, 000 will be irrigated by water diverted at the Timba damsite. The un- | 


a regulated flow is sufficient at all times to meet irrigation water requirements. a 
Taking the power output corresponding to the rule-curve operation as a 
| norm, the optimum combination of Timba reservoir storage and the channel Pil. 7 : 
improvement will at giving the lowest cost of of 
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Average annual generation of 346,000,000 kwh 


‘Table 3 presents cost estimat 


a fo: 


reservoir storage. For case a low for irriga- 


7 tion would be constructed at an estimated cost of 22,000,000 pesos. Be-- 


cause of t the low head available and poor streamflow regulation, it would 


of 60,000 kw capacity would wihdas to be “constructed with costs as indicated 


7 Sa . Timba reservoir operated tor maximum power output with flood con- 

trol incidental. _ The firm capacity being 67,000 kw or 7,000 kw higher than © 


> the norm, a credit of Ps 11 million is given in line 11 representing the cost 


within the 85,000-ha area subject to direct flooding fromthe 


_ of 7,000 kw of thermal capacity; in this case, the Ps 11 million actually _ 
represents ¢ cost savings at the already - ~existing ‘Yumbo thermal rmal station, 

3 reservoir according to the rule curve in Fig. 5. “This 


: 


4, - Timba reservoir operated for maximum flood-control effect with 0 
= power incidental. With the firm hydro capacity at 43,000 kw, iy, 000 a 
in _ kw of thermal capacity is poguices to meet the norm. "From line 11, 


at La Balsa and Juanchito for four ca cases 


zz _ as given in Table 4. Figure 7 presents the results of detailed computations of y 
the effect of the Timba reservoir on the 1949-50 flood (Fig. 6) | if mayest dace 


itl The capital cost of the channel improvements required for each 0 of tee t four 
a cases, with design capacities at La Balsa and at Cali as given, respectively, 
a the second and fourth columns of Table 4, are given in line 8 of Table 3. = 
_ _ The total capital cost (lines 12 or 23, ‘Table 3) is increased in each case 
“by the capitalized value of the annual cost of operation and maintenance, _ 
depreciation and fuel consumption. ‘The overall cost thus obtained provides 
the basis for comparison. " It is the lowest in Case 3, which vecpiaa repre- 


high the wet season for a 


- Cauca River and in an additional area of about equal size immediately ad- 
_ jacent to the flood zone on each side; flood control of the Cauca River will 
thus benefit directly about 170,000 ha or over 40% of the valley. . Drainage ae ce 
at works, principally canals and pumping stations, ' will of course be necessary Soy 
to supplement the main flood-control works (Timba reservoir andCauca tt 
‘Se Flood control and supplementary drainage works within this 170, 000. -ha ‘a Be 
al low zone will also benefit the remaining 230,000 ha in the valley but in an in- 
sal direct way. Irrigation of these 230,000 ha is proceeding on a gradual basis. 
_ (The rate will be accelerated « once the Timba dam and the main irrigation — ioe 
a leading from it are completed. ) A deterrent to more rapid develop- Peni. 
ment of irrig ion has been | lack of 0 outlets and 1 high 
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FOUR ALTERNATIVE ‘SOLUTIONS 


Case 1. Case Case 3. 
used for used for 


Timba hydre f 
60,000 


60, 000 


Thermal 
rae Total 


9. Diversion dam for irrigation 
10. Timba multipurpose project 


Thermal station 


ANNUAL COSTS ‘(thousand pesos) 


0 

150 


17. Total 

otal 2,560 

19, Timba multipurpose project aa 


COSTS (million pesos 
23. Capital costs (line 
24. 0 & M capitajized (g) 
Fuel capitalized (g) 
Total 


| 
han 
cost 
ly, q 4 
N 
se 
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Includes 73 million surcharge storages 


down to El 100 
(c) eservoir drawn down to El 1000. 


Include interest during 
to pesos at rate of 1 US dollar = 7 pesos. 
| Depreciation assumed negligible for the channel improvements. — 


“Depreciation computed on sinki ng-fund basis interest 


at 4%. Resulting annual rates are: 


‘ » 
multipurpose 


Thermal station | 40% baie 
coal price 30 pesos per 
metric results in 1.5 centavos” per kwh 


CVC can borrow local rate of 8% and 
funds at 54%. Operation and maintenance, and fuel 
osts are in local currency and are therefore ‘capitalized 
- (8%. As 2/3 of the cost of the diversion dam and of the 
Timba multipurpose project. are in foreign currency the 
average interest rate is 6.33%. For thermal 
stations: 4/5 of the cost ; is in foreign currency, giving a La 
- weighted average interest rate of 6.20% As these rates 
are close to each Baia — a uniform rate of i 


roi 


- 


— 
— be 
3 
— 
— 
— 
wo 


9 aunos 


in 
q 


id 


ot 


al 
< 
~ 
< 


Fang 


1 


ber, 


ptem 


JIOAJESOY 


pik 


ve, 


i | | || <= | 
— \ | \ | 


No channel channel 


2 4,200 nas 190 210 
4 


670 


Reservoir operated according» to o rule curve in Fig. Fer re reservoi 
me levels at or above the rule curve » release from reservoir is such 


» 


CONTROLLED 


Controlled discharges for average ‘flood are for Case 3 
than for Case 2 because of rule-curve operation whereby reservoir 
ae is evacuated in antici ation of possible occurrence of the design ; 
flood (that of 1949-50); but note thet Case 3 provides NEN: 


; po in adjacent downslope : areas. In many cases, a rather anarchic : situa- 


@ has developed in which upslope landowners have thrown their irrigation 


wastes on their downslope neighbors; or, drainage - -ways have been con- 


1949 


4 


— was stated above that drainage works alone would materially 


her damage caused by Cauca River floods. The Cauca River rarely remains _ = 


above flood | stage more than two weeks (many crops, if sufficiently mature, E - 
; can withstand flooding for this long) but the low areas behind the riverbanks 
, are now flooded for considerably longer periods, owing to the absence of ie, | 
4 drainage outlets. ts. However, even if drainage outlets are constructed : as a first ae 
a stage (providing sizable benefits), general flood control will nevertheless _ 
i - goon become necessary to protect the high- -value crops that will be grown, as 
as the drainage facilities themselves. Moreover, most farmers will not 4 


FIGURE 


risk planting of valuable crops in the absence of a sufficient degree of control 
= _ With respect to the tributaries | the situation is different because the floods 
. F are of short duration with little volume in the peak. , Studies of several _— 
- cal areas have shown that the most economical design is to provide a rela- 
a _ tively low degree of protection against the tributary floods with occasional — 
overflows: (once in 4 or 5 years) removed bya system 1 of drainage canals. 
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The Timba Project 


three spillway gates, each 13 m wide by 11.5 m high; a 60,000- kw 
e powerhouse; and miscellaneous facilities including a flood- control outlet alll 
irrigation | intakes. discussion herein will be limited to the spillway 
r determining the. economic justification of 


Spillway D Desig 


4 
7 
the one hand, it is. known that the upper ‘Cauca | basin is ina a region 
q cold polar air masses. On the other hand, the period of actual streamflow a 
records is short even though the qualitative historic record since the six- i 
teenth century. indicates the absence of floods of the order of magnitude re-— 
corded in the United States. The hydrologic problem is further ccingticates — 
_ that it is practically impossible, owing to the paucity of data an and the spotty | 
: and erratic character of storm rainfall, to apply methods now having standard» i 
acceptance involving “storm transposition” and the “unit 
methods were therefore described as follows: 


1. Correlation between basin rainfall and daily discharge were 


4 attempted. The basin rainfall was determined by the Thiessen method; a : 
more accurate method involving the plotting of isohyetal maps was not 7 
7 possible owing to lack of sufficient data. A correlation between the aver- ; 
age daily flood flow at La Balsa and the daily basin rainfall was found to 
ree be unsatisfactory; _ however, a better correlation was found between th the 
total basin rainfall for the preceding five e consecutive days and the aver- a 
age daily flood flow. Considering all recorded daily flows higher than 300 i 
— m3/s, the correlation plotted as shown in Fig. 8. It is thought that the c- 
=e curve » represents the maximum daily discharge that can be 1 reason- : 


ably expected as a function of the 5-day rainfall. 


2. The 1000-year 5-day precipitation was obtained by extrapolating the 
18-year rainfall record up to a frequency of 1,000 years giving a 5-day i 4 
basin precipitation of 186 mm (7.3 inches), or 62% more than the 10-year 

’ storm (see Table 2). It is believed that the upper Cauca basin, because “a 
an its location in the equatorial latitudes, can receive a probable maximum 
=, depth of rainfall of no greater order of magnitude than the 1000-year rain- ace 


fall. Such a premise is considered sound based on the following 


= 


Smeteorological accidents” in which, however, one or more of the follow- 
conditions are always involved: a tropical hurricane, a . meeting be- — 
tween a cold polar air mass and a warm tropical air mass, or snow- romnelt. 


: “which indicate the expected rate of rise and fall of floods, a hydrograph_ 


was sketched based on the maximum daily average discharge determined 
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resulting hydrograph had the following 


flow 400 m3/s to 480 m3/s 


Volume (in 3 Sdays) 685, 000,000 m3 125 mm 


Comperionn, of the spillway design flood with past floods of record indi- we i 


cates that with respect to peak, the maximum flood of record has a peak - 
dl = charge of 1,225 m3/s; thus the spillway design flood is almost five times as - 
great. with respect to volume, the is made: 


Maximum 1- “day volume 
ast The peak discharge | of 6,000 m 3/s corresponds to a Myers Lb c of “" 
= 550 in the formula Q =C VA a2 where Q is the discharge in cubic feet per - . 
second, and A is the drainage areain square miles. j 


al With the inflow flood established and with the top of dam set at El 1029. 0 — 
4 (based on i ee termes mad the capacity and dimensions ofthe 


Outlet functioning of the t two 

Study of winds (mild in the \ upper Cauca basin) and reservoir - fetch indi- 


cated that a freeboard of 1.5 mis ample. However, owing to the fact that _ bein Fa sf 


- Timba is an earth dam and that its safety depends on proper operation of a 
i _ spillway and outlet gates, it was decided to allow an extra meter of . <<a. 
- ‘board in Condition 1. In Condition 2, it was assumed that one of the two 
operating gates controlling the outlet conduit would be stuck in a closed ee = 
tion; however, as this is only a remote possibility (the gates will be operated eo 
: frequently) a reduced freeboard of 1.5 m was accepted for this case. La 
- Conservatism is also involved in the adoption of an initial reservoir eleva ; 
tion of 1025.0, which is the level of the top of the spillway crest gates. The a . 
_ record of greatest floods, Table 5, indicates that floods of appreciable volume 
: (more than 10% of the flood storage of 305 million m3) have never — 
each other within a short time; nor have such floods occurred from enrly iJ 


adhered to, , the spillway design flood is not likely to occur on top of a full 


reservoir. ” On the other hand, it was felt that pressure from farmers along Syl 
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CAUCA VALLEY 


TABLE 

PLOOD VOLUMES ABOVE RESERVOIR RELEASE 


Streamflow records, October 1945, to 


aio? Volume Above 


‘Mey 3, , 1949 20.1 


1950 


a 


+a 
“1954 Bes: 931 
Jens 1) 105 1956. 1,225, 
aw Reservoir release | of 650 3/5 would be about 
average maintained ‘while storing the 
_ the river could, in the future, upset the operating rules; they might insist - 
s that reservoir releases be kept lower than the amounts. planned, even though 2 


The criteria in a of three radial gates 


ae The two. principal elements of a balanced flood- ‘ebaitrol plan are the Timba z | 
‘reservoir combined with the Cauca River improvement (Case 3 described — * 4 


above). . Planning and design are well advanced for the Timba reservoir but 

not for the Cauca River improvement. Preliminary design of the latter has, — 

however, been completed as summarized in Table 6; in general, design 

capacities are 20% to 25% higher than existing channel capacities. © a 4 
q tedly approximate, the preliminary design is still considered good enough for 5 4 

the making of economic comparisons, which are on estimated annual 


“costs of alternative long- range plans. 


river improvements) will be a highly 
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CAUCA RIVER CHANNEL IMPROVEMENTS 


Present River 179 
Future River Length 
River Shortening Due 


Length _of Right-bank 


“The 
‘meanders are most ae ia in the | upper 100 km where the prin is 
- greatest; this is apparently explained by a higher proportion of bedicad 
and higher erodibility of the riverbanks in this portion. ‘The riverbanks, 
=< 
in general, consist of fine - -grained silty sand having relatively low an 
aa to erosion; the river bed consists of sand and gxavel. hae : 


Since ‘records began at Cali (Juanchito Bridge) in 1945, the river has 


been degrading at an average rate of 5.8 cm per year. on ‘(This degradation ‘ 
os be due, in part, to dredging of the river for sand and gravel to supply a 


the local construction industry which has been increasingly active in re- 


ac. a cent years. It is estimated that 350,000 m3 are removed annually by hand, 


a about 100,000 m3 by mechanical means, withinareachofriver | 
several hundred meters long just upstream of of the Juanchito bridge >near 
Cali.) Such degradation has not, however, | been noticeable at points along 


_ 3. The Timba reservoir will have the effect of removing a major ae 


ee of the sediment (suspended as well as bedload) now transported by ae 
=. river, thereby increasing its erosive power; on the other hand, the re- eae 


ie ‘Sediments will continue to be « carried to the main river by the tribu- 


‘taries.. i However, some of the largest tributaries, , such as the Palo, Fraile, 
Bolo and Amaime, traverse the 1 valley floor for a distance of about '25 km; 


typically, these tributaries overflow and deposit the bulk of their ieee | a od 


— 
| 
9, the Cauca isan 
— 
— 
— 
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i CAUCA VALLEY 
alluvial on the valley floor before reaching the river. 


» reachin 
_ 5. The tributaries vary in their sediment-carrying characteristics ow- 


~ ing largely to differing soil characteristics and vegetal cover. . For faite 
a the Rio Bugala-grande is a notorious silt carrier; a —< fh 
- conservation program within its watershed would produce benefits to local - 
"agriculture and, incidentally, would greatly reduce the quantity of sediment 
reaching the Cauca River from this source. 


. At present, the river overflows frequently and deposits 
quantities of sediment on the flood plain. In t the future, after construction 
of Timba and the channel improvement, the net sediment- load remaining y - 


iu will be confined to the channel instead of being deposited on the flood plain. q 
7. Be 


low Cartago, the river changes from an alluvial to an ‘erodible- bed 
os, flowing for the greatest part of its length through relatively nar- — 
- _ row canyons to its mouth at the Magdalena River. There is practically wd a 
settlement along the river below Cartago; therefore, any of the 
wah) 8. Excessive degrading of the river in the reach from ad to Yumbo al 
would create serious dislocations because of the existence of several im- 
_ portant water intakes of which the most important are those for the munici- 
pal water supply of Cali, and for several important ind industries | near ‘Yumbo — 


including a large steam-electric station. "ee? a 
9. A vast hydroelectric project (capacity abo about 4 500, 000 kw),’ for wh which ich 
a reconnaissance study has been completed, involves diverting about 1: 130 ‘of 
m /s continuously from the Cauca River to the Pacific slope by means of 
‘y a tunnel through the Cordillera Occidental. This diversion would require 
a barrier dam near Vijes. ‘Such could 


Some channel degradation in other aidiihe of the river or would & be 


— as there would result savings in the cost of levee construction | = i 


and improved drainage conditions for land both within and adjacent to the 


. For excavation of the the “pilot cat” method is planned 
, aa most notably on the lower Mississippi where the original river length P 
_ of 680 miles from Memphis, Tennessee to Baton Rouge, Louisiana, was co. 

E - shortened 25% through use of cutoffs). _ Techniques to be followed in apply-— 


ing this method have been described by Matthes. (2a, 


12. The estimated average height of the levees (Table | 6) is 2n m 


3 4s La Balsa to Buga and 1.5 m from Buga to Cartago. “Levees will generally — 
i be set back from the river-bank a distance of 25 to 50 m in order to pro- S 
: vide a margin of safety against bank erosion and to carry out erosion- let 
ae — measures within the 25 to 50-m belt, to consist eeeeiiedind of vari- 

ous types of with root 


a 
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ibu- 
_ As described in the intr flood control will be 
that resulting from more ir rather than from prevention 


of to existing crops improvements; this benefit should be 
ured as the difference in its productive capacity under conditions expected 8 % 
__ with and without flood control. (4) _ To estimate the difference in productive % 
capacity it would be necessary to determine the net income expected from the 
¥ land and from associated resources with and without the project. However, 
_ because of time limitations and lack of specialized personnel, an alternative - 
_-method was used based on estimates of the increased market value of the af- _ 
fected land areas. Checks on these market values were obtained by analyzing 
the productivity of the land for typical crops with and without flood control. 
_ Asa first step, net farm income per hectare for various land classes wa a 


; estimated in Table 7 on the basis | of rental value plus added value produced > 
“ie by management. Rental values were taken at 10% of the selling prices for 
land in the Cauca Valley, for which conservative values, somewhat lower 7 
than currently prevailing prices, were used. (The rate ‘of 10% is the current 
interest rate in Colombia for relatively large private enterprises. 


_ Added value produced by management was taken at 30% of rental 


Entrepreneurs customarily pay landowners a rental equal to a third of the 

_ gross profit. The remaining two-thirds is for entrepreneur’s gross profit — 

_ including bank-interest charges, management services, income taxes and net 

profit. The latter two items are conservatively estimated at 10% and 5%, 

respectively, of the entrepreneur’s gross profit; this is equivalent to 20% 
_ and 10% respectively of the rental value; see note 2, Table 7, Income taxes, 

while not a part of net farm income from the entrepreneur’s point of view, — 


: wealth produced by the project. and should therefore be included in _ 


_— check on the ) assumed selling price for Class G land is provided by con- 
an oe current values and d production costs for rice and soybeams, typical 


_ es i The present average price of unmilled rice per sack of 150 pounds 
(avoirdupois) is Ps 55; for Soybeans, the of 125 5 pounds is 


4 
2. Ina 3-ye -year r cycle, at least four harvests a rice and two of beans a 
- obtained. The yields per plaza (0.64 ha) are 47.5 sacks per harvest of : 7 
Tice and 16 sacks per harvest of of beans. T The resulting gross | 
. Production costs per per plaza are Ps 620 Ps 
_ 250 for bea beans, which gives an average production cost of Ps 1,000 a 
plaza annually. All production costs are included except entrepreneur’ 
interest charges, the value of his management services and 
a. these are considered part of entrepreneur’ Ss gross profit (see item 


5). "Gross profit i is therefore Ps 3,000 per plaza. 


Allowing 10% for a area by roads and canals, 


The remainin Ps 1,800 is entre reneur’s ross rofit. PSH 
g ntrepr ‘gross profit.) 


ae a 6. As use of ‘money in n Colombia now costs about 10% per annum, Class 2 


G land is worth Ps 9,000 per plaza. To be aE a value of Ps a 2 
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ANNUAL FARM INCOME PER PER HECTARE 
(s 
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1. Rental value taken at 10% of sale pr ise. 
t of Added Value Produced by Management tine at 30% of rental value 


< Net Farm Income "Bt equals Net Farm Income "A" less annual cost ef irri- 
_ gation and drainage works, taken at Ps 150 per plaza (applies only te ; 
_ Glasses F and G). The latter cost is (i) 10% of the estimated con- 
struction cost of Ps 850 per plaza of these works limited, however, te 
canals and drainage ditches (main canals and pumping stations 


m 


* 

a 

if 
con- 
ids 
a 
— 
fe. 
oon (a) Reduced 10% to all 
the (b) One plaza equals 0 
ogg 


The latter area is subdivided into land ‘classes, with | 
ing benefits per hectare, in Table 8. _ Also indicated in this table are the bene- : 
fits in the areas adjacent to the flood zone due to lowered groundwater condi- a 

tions. The ultimate direct flood-control benefit, in Part I of Table 8, is Ps 26 E's 
million from La Balsa to Buga and Ps 12 million from Buga to Cartago or = i 
total of about Ps 38 million. _ Following construction of Timba and the river _ 


impr ovement, the completion of the bulk of which will take 4 years, full culti- a 


=» 


vation of the benefited lands will require additional time: this has been as- - 
gumed to be 5 y years. The ultimate direct flood-control benefit has therefore 
been discounted from Ps 38 million to Ps eae a as shown in Part Il of 
Additional indirect flood-control will value added through 
processing of increased agricultural production; reduction in loss of human | 
life; improved sanitary conditions; reduction of loss of livestock by drowning; — 
= savings in cost of highway and airport ‘construction and maintenance — 
from improved flood and drainage conditions; and elimination of flood Gamage 


oie _ caused by interruptions to « commerce and to industrial production. Although — 


ig om least half of the direct benefit. The total flood-control benefit then is 1.5 
‘The annual cost of flood control was determined by adding the capital cost 
4 of the river improvement to that part of the Timba project capital cost al- 7 
located to flood control, as estimated on the basis of the “justifiable- -= api 
 alternative- ~expenditure method.” It was assumed that the foreign- 
part of the resulting total capital cost will be amortized in 20 years and that Fe 


_ then made for depreciation | and d for operation and maintenance ot The resulting 


interest on the local-currency part will be paid indefinitely. Allowances were : 


benefit of Ps 48 1 is about 4.3 times the annual cost of 
Ps 11. million, the two projects are clearly justified. 
eS - For planning purposes it has been assumed that complete development _ 
~ (including irrigation and drainage, as well as flood control) of the entire val- 
ley floor of 395,000 ha will be accomplished during a 30-year period. The © 
_ key structure in the program is the Timba dam. Owing to its cost and multi 7 
_ purpose nature, several years are required to arrange for its financing. ey z 


. However, in the meantime, it is possible to proceed with the reclamation of ae 


dependent on the Timba reservoir for irrigation | water. Moreover, these 
i initial projects may be regarded as “pilot projects” as they will point the way a 
— efficient development of the bulk of the area in later years. ae 2 
Three such pilot projects are Roldanillo- La Unién-Toro he 200 ha), 


Aguablanca 


ion 
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Benefit 


3,2 


2 __700 20 1000 
45 19,800 35,000 9,900 


in. 


‘Total Uieinate Benefit Ps 37, 


(a) Difference in net farm income between Class 


‘Estinated benefit due to lowered groundwater in area acent 
J 


Assume 5 years will be required to. sshheve the ultimate direct benefit. 
interest rate of 8%. Determine re value of the benefit for the 50 
rs following project completion. 


) Present value of an Annuity © Ce rtain for 45 years i is 12. ny; 


2. ll x .681 = 8.247 
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er, 1950 

_ Was | based ‘mainly on: extent of local interest, technical feasibility and econo-— : 


mic return. _ The first two projects : are under construction while the third is he 


the land reclaimed by this project ‘will ultimately t be to 


a Ss oe (see Fig. 9), the degree of ‘protection will be much higher than that adopted | 
_ os _ elsewhere | in the valley. ‘The river levee on the east is being built to a profile — 
2.5 m higher than the highest recorded—during the flood of February 1950. Re 

The ‘frequency « of a 1. ‘5 m higher (allowing 1. m freeboard) is estimated 


a a Elimination of flooding within the e Aguablanca project area will eliminate 55 - 
million. m3 of natural valley storage or 7.5% of the 740 million m3 between __ oe 
ne _ La Balsa and Cali. _ The peak discharge of the 1950 flood was 1,000 m3/s. The ‘ 
_ eliminated valley storage would cause an increase to 1,032 m3/s, for which m 
% the river level would be only 15 cm higher than under present conditions. a 
_ Although, as described above, diking of the entire flood zone in the valley 
is not feasible, it was decided nevertheless to make an n exception in the case 
of Aguablanca ‘for the following reasons: * (a) The v valley storage eliminated by 
- this one project has a negligible effect on river levels. (b) Benefits will be 
_ much higher in the case of the Aguablanca project, per hectare of reclaimed — 
~ land, than elsewhere in the valley. (c) Further urban encroachments on = 


flood plain are not visualized, at least for the next generation or two. (The 

_ remaining towns are all located well outside the flood zone. oh . eee 
The interior drainage system includes (1) a pumping station to remove in- _ 
erior drainage water during flood stages in the river and to lower the pres- — 
_ ently excessively high groundwater level; and (2) a dredged storage pond to 


reduce pumping requirements and, at the same a 


Design data with to the pumping station a are as follows: 


area 


of drainage area 


years 
inflow 
Peak outflow 


-year flood (Design Flood) 
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oe terior area with a 3-year flood in the Cauca River. “The » available data intl 
- = cate that these two events, if not entirely independent, are almost so and that 
the of ‘coincidence is about once in 7 years. 
For design conditions 1 sitet river r level at 953. 8, static heads will be 4.3 “- 
maximum and 1.9 m minimum. . However, pumping will be required also for 
river levels equal to the sump level or as high as } 956. 0. The minimum static 
head is thus zero andthe maximumis6.5m | 
-_ _ Because of the largely independent nature of the two hydrologic : factors z 


sirable, asa second ‘requirement, to have a pump with a curve 
of head vs. ‘discharge in approximately inverse relation. 
To best meet these two requirements, vertical propellor pumps with vari- 
_ able blades are being specified. This type has the further advantage of low a ; 
starting current, which is important at the Roldanillo-La Union-Toro project 
_ where pumps identical to that at Aqui will be used and where — 7 
Although ‘unavailable in the United States in large 
sizes and on special order, variable-blade pumps in relatively small pa 
27,000 - -gpm pumps will be about 42 inches in « diameter) are now being 
“¥ ‘manufactured in the ‘United States and Europe by several well-known firms. , 
me Bids recently received have a that the variable- blade type provides 
_ the most economical solution. 
pondage of 770, 000 m3 will composed of 600, 000 000 m3 in ina storage 
- pond and 170,000 m3 in the drainage canals. . If. the drainage canals are later 
_ replaced by conduits or lined channels, thus increasing velocities and reduc- 
storage-in-transit, the storage pond 
oa oa 2 To obtain the 600,000 m3 of usable storage in the storage pond, 840, 000 mo 


1 Without the storage pond, of the pumping station 4 


"pond. If the capitalized cost of depreciation and electric energy ‘consump 7 


the economic comparison is in favor of the ‘storage pond. 


the interior drainage area becomes over the course 


‘oO? - the years with storm sewers and streets constructed, runoff characteris- 4 
ee ties will become progressively worse in the sense that the volume and es- } 
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ncrease four is, up to 36 to 54 m3 as compared with 
18 m 3/ under initial conditions. the storage pond, little if any in- 
a in the pumping capacity will be required to meet the needs of in- _ 
 — urbanization. However, to be safe, space is being provided at the 
pumping station for two additional pumping units, when and if needed, which — 
would raise the rated d capacity from 81, 000 to to 136, 000 gpm. 


aan _ The value of the general area in which t the par will be located will oa q 
be greatly enhanced. This area is now a swamp, 62 hain extent. The 
Pe eZ pond will occupy only 46 ha of which the water area will be 31. = 
ha and 15 ha will be in 1 islands. Moreover, the area ‘surrounding the stor- a 
_ age pond is being reserved for parks and housing developments. The ex- — 
- cavation is planned so that the : shores of the pond will be well defined, = 


=e flat shallow places. A minimum depth of 80 cm is being provided ep > 


law adopted by the in April 11958 provides thai 


: cvc in n the form ¢ of an 1 assessment which c can 1 be paid within optional periods ss 
of either 5 or 10 years. For the shorter term, interest of 8% will be charged 
7 and, for the latter, 10% al (The u usual interest | rate in Colombia for commercial 
4 credits is 10%.) In addition, landowners must pay CVC 10%of the enhanced a 
value of the land benefited by the project, as determined by a committee of | 
three assessors chosen jointly by CVC and the landowners. 
a Although many landowners are opposed to this formula, 1, they should have no — 
‘difficulty in paying the assessment it calls for. Unimproved flood-free land we * 
= the project area has a market value of 15,000 to as high as Ps 50,000 ie ; 


per ha with an average of about 20, 000. In value 
tains a provision whereby may elect to pay of their 1958 con- 


project benefits, as measured by the enhancement value 


Ps 15,060,000 giving a benefit-cost ratio of 4. 1. 


Roldanillo- La | Union-Toro Project & 


estimated project cost including -of-way is 


An of this project, Fig. is the nature 
2 two of its features: the La Union Pumping Station and the intercepting oo" 
| ie The La ‘Union Pumping Station will contain three 27 ,000-gpm units to pump 
_ drainage water to the Cauca River during flood stages; also, to pump water 7 
: from the principal drain to irrigation canals located on the high bank of the 4 ‘ 
river. As static lifts are similar to that at the Aguablanca project, it og 
been possible to adopt identical me the at at both» 


_ handle from s streams draining the east ape of the Cordillera 
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western of the area. accomplish this second purpose, the 
will be provided with check structures each including a tainter gate operated | 
"automatically to maintain water levels sufficient for irrigation; during floods, 
the gates will open completely so as to provide unobstructed passages. 
‘The natural levee along the left bank of the Cauca River provides protec- se 
~ tion against floods up to an estimated frequency of once in 5 years. When the 
_ irrigation canals on the Cauca River bank are constructed, the excavated 
_ material will be placed to form a levee 1 to 2 m high; the protection then 7 
_ available will be against about a 7-year flood. This degree of protection is — 
0) 
F considered adequate in view of the low degree of present land use and intense ss - 
desires of the small landowners in the area to cultivate their land. _ The future | 
Cauca River improvement will increase the degree of protection to about 10 
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FRC TION FACTORS IN CORRUGATED METAL PIPE 


more 


study made to friction coefficients in large ‘corrugated 
metal pipe by measurement of head losses in pipes 3, 5, and 7 ft in diameter a 
and to determine velocity distribution f for ' several discharges in each test — 
so that the results could be generalized to show the relationship between fric- 
tion coefficients and absolute roughness is discussed in this paper. Data taken oy 


during the study ar are presented tc together with recommended values of friction -_ 
and roughness coefficients. 


factors for various types and sizes of pipes and conduits. frovsecincey compara- _ 
) tively few tests had been made prior to 1950 on corrugated metal pipe of the Pat : 

larger sizes commonly used in engineering practice. For this reason, — 5 
Office, Chief of Engineers, U. S. Army, authorized tests on 3-, 5-, and 7-ft- SS: 7 
diam corrugated metal pipe to be conducted at the Bonneville Hydraulic | oy | 

‘In 1926, Yarnell, Nagler, and Woodward published results of hydraulic 

tests(1) at the State University of Iowa on 12-, 18-, 24-, and 30- in.-diam cor- 
rugated metal pipe that varied in length from 24 to 36 ft. The results of these 
tests are applicable to similar installations of short pipes and culverts, such 

as those under highway and railroad fills. However, they are not directly ap- 
plicable to long pipelines because entrance, transition, and exit losses cannot a 


Note: Discussion open until February 1, 1960, To extend the « closing date one month, “os 5 
a written request must be filed with the Executive Secretary, ASCE. Paper 2148 is 
part of the copyrighted Journal of the Hydraulics Division, Proceedings of the > ae a-, 
American Society of Civil Engineers, ' Vol. 85, No. HY 9, 1959 
1. Chf. Hydr. Design and Lab. ‘Section, U. Ss. |. Army Engr. Dist., Portland, 
Corps of Engrs., Portland, Ore. 
2. Hydr. Engr., U.S. ‘s. Army E Engr. . Dist., Portland, Corps o of Engrs., 
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length in 1 relation to their - diameter. Keulegan(2) used velocity 


_ data from the foregoing tests to compute friction factors for long pipelines | 
for pipes of larger diameter. = 

In 1950, Straub and Morris published results on . studies(3) conducted at the | 
St. Anthony Falls Hydraulic Laboratory, University of Minnesota, of flow __ 
_ through» 18-, 24-, and 36- ~in. -diam circular and arched corrugated metal pipe 
193 ft in length and laid on a slope of 0.002. These studies showed a definite 
a systematic increase in the values of the friction and roughness coefficients _ 


= and “n” with increased Reynolds numbers | (Re) over the entire range of 


- most ‘ “f” vs Re curves for both smooth and rough commercial pipe have 
shown constant or slightly characteristics throughout the experi- 
- mental range. . Values of “n” from 0.0216 to 0.0252 were obtained for cor- * 
rugated metal pipe flowing “full and from 0.0228 to 0.0258 for open-channel 
2 - Values recommended f for use in design were 0.025 and 0.0235 for full- 
. The purpose of the tests at the Bonneville Hydraulic Laboratory was to "e 
on friction factors for 3-, 5-, and 7-ft-diam corrugated metal pipe as 
indicated by head-loss measurements for a range of velocities } up to approxi- 
— 10 fps for 5- and 7-ft-diam pipe and 16 fps for 3-ft-diam pipe. Tests 
“= made with both full-pipe and open-channel flow in the 3- and 5-ft- -diam = | 
pipe and with full-pipe flow in the dae ft-diam pipe. The initial test schedule Hy, 
was ‘expanded to include studies on the 5- and 7-ft-diam | pipe with 25 per cent 
of the periphery (at the invert) smooth paved and on the 5-ft- diam pipe with — 
50 per cent of the ‘periphery smooth paved, 


Description of Project 
aa _ The tests were conducted at a site located on the Washington shore. of the 
: c Columbia River at the north end of Bonneville Dam. _ This site was selected — 
i because sufficient space was available for the test and metering sections and ~ 
a because | sufficient head and water supply were available to produce the dis- mel 
_ charges and velocities required for the tests. The layout for the 7-ft-diam 


pipe, typical 1 for | the three sizes of pipe tested, is shown in Fig. 1. The 3-ft- 

_ diam pipe in the test section is shown in Fig. 2. . Water for the study was 
taken from one of the fingerling by-pass channels through the main dam 
: From the fingerling psa ring the water flowed through a 6-ft- -diam, smooth- 


a 7-ft-diam, smooth-wall pipe and then dropped 11 ft to the metering section 
where it was measured by a system of flat-plate orifices and then returned — 
— to the river through gate valves and short outlets. The 11-ft difference in od 
i elevation between the test and metering sections was sufficient to Re 4 
adequate head on the meters for the tests with open-channel flow. ed = 
corrugated metal test pipes were of standard commercial manufac- 
dle ture and arrived at the test site in 20-ft-long pieces composed of a series of 
: 2-ft-long overlapping sections. The sections were riveted together and each 
contained lapped and riveted longitudinal joints that were staggered so that no 
two were in line in adjacent sections. Before installation, the riveted joints 7 
| were scaled against leakage by welding on the outside of the pipe. The 20-ft 
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grade and alignment along a continuous footing between concrete saddles. _Ob- 
servations were made when the pipes were : flowing full but no sagging or de-— 
= was visible. The average diameters (top of corrugation to top of ; a 
- corrugation) of the 3-, 5-, and 7-ft pipe were 2.996 ft, 4.95 ft and 7.05 ft, — 7 
respectively. depth and longitudinal spacing of the on 
pipes were 1/2 in. and 2-2/3 in. , respectively. 
The invert paving of the 5- and 7-ft-diam test pipe —_ accomplished by 
ine 18-gage, 10-ft- -long sheets of metal over the corrugations 0} on the bot- 
tom of the pipe (see Fig. 3). "i The lower sector of the p pipe was then coated with 
| bituminous hy = complete the simulation of the conventional asphaltic on 


gor 


‘apnea system of flat- -plate orifice meters designed to measure a continu- 

ous range of discharges from 8 to 385 cfs and to provide an overlap between — 

_ meters operated with a differential head of from 1 to 9 ft. An aperture-plate 

_ diameter ratio of 0.7 to 1.0 was selected for all orifices except the one in the 
_ §-ft section; this orifice had a 0.66-to-1.0 ratio to secure slightly different 
characteristics. ‘The orifices were machined from steel plates that 

varied in thickness from 1/2 in. for orifices up to 3 ft in diameter to 1 in. for 
those larger than 5 ft in diameter. Piezometer taps were located one pipe 

diameter upstream from the orifice plate and one half the ‘orifice diameter 
tc downstream from the plane of the orifice. _ The meters we were assembled prior 4 
to rating and were moved intact to the metering section. Periodic inspections 
Were made of the orifice plates and piezometer taps to insure that they re- ‘7 


mained clean and free from rust. In the metering section, the orifice meter 
a the 6-ft-diam pipe (Fig. 1) was preceded by 19 and followed by 6 diameters 
‘- of straight, constant-diameter, smooth-wall pipe; in all of the other meter 
installations, the approach and exit pipes exceeded the length standards pre- 
x scribed by the Special Research Committee on Fluid Meters of the — 
st Society of Mechanical Engineers. (4) First, two 24.0- by 16.2-in. orifice 
meters were calibrated separately, u using a standard calibrated venturi meter, 
E and then installed in the metering section. Successively larger orifice meters 
were calibrated as follows: the entire length of constant - diameter, — 


that contained the orifice were placed in the test section, the ebitine piezo- 
ie meter taps were attached to a standard differential water | manometer, om 
varied rates of flow were introduced into the system. _ The quantity of flow, 
_ as indicated by previously calibrated orifices in the metering section, to- 
roe = with the observed head differential on the test orifice, provided the 
: data from which experimental coefficient curves were plotted. After i 


all data were checked for accuracy and consistency, each orifice meter and 


adjacent — were moved to the metering : section and installed for further 


. ; ot The test sections in the 3- , 5-, and 7-ft-diam pipe were 180, 251, and 80 
long, respectively. Pressure observation points for the T- -diam pipe are 


/ indicated in Fig. 1. There were 12: stations in the 3- ft pipe at approximately 4 
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spacings and 10 stations in the pipe a at ‘approximately 20-ft 

Each pressure observation point consisted of two 1/8-in. -diam piezometers 
that were drilled into the valleys of the | corrugations 45 degrees above the 
invert on each side. The two piezometers were interconnected by a common 

lead to a single manometer in the bank shown in Fig. 4. Individual pressure - 
; ‘fluctuations at each piezometer varied from 0.005 ft at low discharges to a 
: 0. 015 or 0.02 ft at high discharges, however, it was found that consistent re- = 
sults could be obtained from a series of 5, 10, or 20 consecutive visual - - 
es was determined ina adie study y on 3-ft pipe (Fig. 5) that these ia 
» _ piezometers were influenced by the partial recovery of velocity head in the 3 
corrugations and therefore indicated pressures higher | than the true static 
pressure as taken at a point 6 in. within the pipe. However, they were con-— # 
a: sufficiently accurate to determine the differential pressure between XM 
is "successive stations along the pipe axis. The data in Fig. 5 are qualitative in 4 
nature and are included to show that the . uniformity of pressures measured he 
_ along the corrugations depends on exact placement of piezometers. It is be-— 
= that the | nonuniform points shown on some of the grade lines presented — 
in Fig. 7 were caused either by abnormal velocity distribution or by — 
; fections within the pipe that affected localized weequities and pressure — 


Static pressures during distribution studies and free water- 
surface levels for open-channel flow were measured by means of tubes that ' 
a a were the equivalent of the static leg of a conventional pitot tube (Fig. 6). The 4 
static tubes were : checked to determine the effect of local flow acceleration. | 
A At first these tubes contained four piezometers that were connected to indi- 
ie vidual water manometers, and the instrument was aligned with the flow by a 
.- rotating the tube until the pressure, as reflected by the manometers, was ting — 
equal at all piezometers. Prior to the measurement of grade lines 
= _ in the 5-ft-diam paved-invert pipe and the 7-ft-diam pipe, the static tubes . 
Were revised so that the four piezometers in each tube were connected to a — 
- common manometer that reacted much faster * pressure changes than did ae 
the instrument oforiginal design, 
- ye Two types of impact tubes were used to measure local velocity heads from 
_ which velocity distribution curves were determined. The type used first re- 3 
_ sembled a pitot tube but lacked the static openings on the sides of the barrel. 
_ The tube was mounted in a bracket that slid along a strut placed across the b 


| ; } pipe diameter. The impact tube was satisfactory during tests on the 3-ft pipe, — 


is although at high test velocities it vibrated and was difficult to keep pointed — 4 ed 
into the flow. It was found that the tube would require additional support © : 


brackets when velocity heads were measured in the lower half (extended posi- a 


_ tions) of the 5-ft pipe. Because it would be necessary to unwater the pipeline — 


‘ ee during each run in order to install the additional brackets, the impact c cylinder = 
(shown in Fig. 6) was substituted ae the impact tube. This was sup- 


_ manner: after the flow was ee to a desired rate through « one or more of 
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meters and valves, and the average differential 
heads on the meters were derived from a series of 10 readings, the two ob- — 

_ servers proceeded to the test section and measured the hydraulic grade vse 


within the test pipe by means of 10 consecutive readings of the water- surface 

_ elevations within the manometer tubes. Then the observers returned to the -. 
_ metering section and checked the meter ’ readings. In all tests, if the average 


- differential head on the meters at the beginning and end of a test run agreed an 


of flow; otherwise the entire run was discarded. 
_ _ Velocity distribution in the 3-ft-diam pipe was measured by means of a 
i combination of one of the previously described static tubes and two impact _ ‘” 
_ tubes at each observation station. Leads from the three tubes were connected ag 
- to a 3-leg differential manometer where velocity heads were recorded a i: . 
' groups of 10 simultaneous readings. _ The static tube was mounted 6 in. within i. 
| @§ the pipe | and 45 degrees above the invert, whereas the two impact tubes were 1 
- traversed through the horizontal and vertical diameters on separate struts. . 7 
h order that the proximity of the horizontal strut might not affect the ac- i" 
curacy of observations, the | horizontal strut was removed while the vertical pay 
traverses were being made. _ The velocity distribution in the standard 5-ft- ee 
' diam pipe was measured by means of three impact tubes and one static tube 
that were used simultaneously at each station. One impact tube was used on e 
each half of the horizontal traverse and the other on the vertical traverse. +S 
Instrumental data were observed 20 times in succession to determine the 
average velocity for each position of the tubes. Two impact cylinders anda 
= tube were used to measure velocities in the S- ft-diam pipe with a paved 
invert and inthe 7-ft-diam pipe. 
‘For tests of open- ~channel flow in the 3- and 5- pipes, the approxi-— 
=z rate of flow was regulated at the intake gate; : and the depth in the test at 


pipe was balanced by adjusting the valves in the metering section until the a” 
_ hydraulic grade line, as recorded by the same piezometers or static tubes nf 
P used for the studies of full-pipe flow, was nearly parallel to the pipe invert — iy 
“heart _ Then the depth of flow was measured by means of a static tube located a 

hear the longitudinal center of the test section and the rate of flow was deter- — 

"mined in the metering section. It was difficult to maintain a constant slope of t = 
“the water surface in the test pipes, and continuous fluctuations of the forebay if 
level made it impossible to conduct tests of open- ~channel flow during certain of 
periods when Bonneville Dam was used to absorb peaks in the regional power — 
demand. Tranquil flow was required to avoid the occurrence of standing wave “a 
in the test section; therefore, it was necessary to decrease the slope of the a 
5-ft- -diam, paved-invert pipe ‘from 0.005 to 0.002. The same change and pos- 
sibly revisions of the intake structure would have been — to permit 
tests of of flow in the 7-ft-diam Pipe. 


ty 


BAe average 1 velocity of flow for b both standard | and paved pipes flowing full 


‘diameters of 3. 00, 4. 95, and 7.05 ft. In the case of the paved pipe, the area 

_ reductions effected by the paving are reflected in the values of “f” and “n.”__ ‘tt 

The values then determined from straight- ‘Portions of the 
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_ a . will be noted that the grade lines for all tests are fairly straight with the | ex- 
— _ Ception of the 3-ft pipe tests shown in Fig. 7. In this case, only the straight _ 
] ; io portions of the grade lines between measuring points 3 and 10 were used to sigh : 
compute the slope. >. The other points, near the ends of the test section, were 
ae apparently affected by entrance and exit conditions. Values of “S” obtained a 
ai from these grade lines were plotted against measured discharges to form the 
a. rating curves shown in Fig. 10 and were used to derive numerical equivalents 
for “f,” “n,” and Re that are shown in Figs. 11 and 12. f rie friction factor dd 
“f” was computed from the Darcy in 


was: from the friction factor by the 


ive Typical curves of distribution within the test pipes are shown in 
“Figs. J 13 and 14. Local velocities: were re computed from the formula Vj = V2gh, 
where was the mean xan difference b between readings obtained from the pre- 
- viously d described impact and static tubes. In the vertical traverse in the 3 62h 
- ftt- Gam | pipe, where the impact tube wa was displaced « downstream from the static a 
4 was applied to the value of “hn” before the local velocity w was computed. Also — 
i so! shown in Figs. 13 and 14 are curves of the relationship V/V, as as computed nt 


from the formula: Vy, = VRSg Vv ‘This relationship was one of the 
means used to analyze t the adequacy and consistency of a er aa 


‘The first ane of “n” values for open-channel flow in 


a 3- and 5-ft-diam pipes was the computation of cross-sectional areas and bd 
hydraulic radii for observed depths and rates of flow in the test pipes. Then — 
ae were determined from straight-line portions of plotted hydraulic grade 


and values of “n” were computed from the Manning formula. The 
_ curves presented in Fig. _ 15 show the variation of the roughness coefficient vb 
with the depth | of flow. _ In order to plot rating curves for open-channel flow, © 
ee observed m re - reduced to equivalent v values fora a — of 
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& proportional t to the re root of the slope > therefore, measured discharges 


by the factor | obtain discharge __ 


quantities for the rating curves. Average flow against 


The tests showed that the value of Darcy’s for standard corrugated 


- metal pipe varied inversely with the pipe size (from 0.077 for 3-ft pipe to 
. 0.054 for 7-ft pipe; Fig. 11) and that the value of Manning’s “n” was nearly 
a constant (0.024; Fig. 12). The values of both coefficients for each pipe size 

ie varied with Reynolds number and reached a maximum within the test range. 


_ The invert paving reduced the resistance to flow and caused lower values for _ 


_ both coefficients (“f” values of 0. 048 and 0.042 for the 5- and 7-ft pipes 25 
7 per cent paved and 0.035 for the 5- -ft pipe 50 per cent paved; “n” values of i, 
0.021 for the 25 per cent paved 5- and 7-ft pipes and 0.018 for the 5-ft pipe 
ae 50 per cent paved). Values of the plotted points shown on the foregoing __ 
; - figures are presented in Tables A to G. Typical curves of velocity distribu- 
_ tion within the respective test pipes are shown in Figs. 13 and 14. — 
_ The value of Manning’s “n” for open-channel flow in standard pipe se | 
_ found to be almost independent of the depth of water over the range tested J 
and was approximately equal to the value (0.024) obtained for full-pipe flow 
(Fig. 15). The values of “n” increased as the depth of flow increased in the 
paved- -invert pipes. ~The curves shown in Fig. 15 were extrapolated down to 
an assumed value of *n” 
Fig. 16 when the flow was contained within the paving. This value is com-— 
used in good the flow of water through welded steel or =| 


ou 


| 
ne over-all results of the ied sihtinenaih the findings of the tests made 


St. Anthony Falls Hydraulic Laboratory(3) that friction factors for a given — 


r ‘size of corrugated metal } pipe e are a function of Reynolds number (Figs. 11 
} and 12) and were in general agreement on values of “f” and “n” within the | 
of that t study. However, the tests at St. Anthony Falls were 
within the range where “f” and * n” “increased with Reynolds number and im-_ 


plied that the coefficients would continue to increase indefinitely, whereas _ 


- this ; Study | determined maximum values of the coefficients for each pipe size. ea 


such as of pipe, it 


a = 0.011 to compute the lower | parts of the curves in 


be remembered that flow in corrugated metal pipe does not neces- 


“ sarily follow the laws established for flow in other types of pipe, —— | 
; _ has established that | “wake-interference” flow is obtained on corrugated sur 
moe 
a faces of commercial conduits, whereas “isolated-roughness” flow occurs ov 
surfaces of concrete, Tiveted wood= “stave, cast-iron, and brick con- 
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for use with corrugated metal pipes flowing full is presented in Fig. 17, in _ 

ie “ which the “n” curve for unpaved pipe was extrapolated to include pipes ‘of ee 
from 1 to 10 tt! in diameter, and values for the “f” curve were computed from. 
values of “n.” Points from the tests at St. Anthony Falls (1.5-, 2-, and 3-ft — 
pipes) are included to show the validity of the curves for pipes of smaller a) 
diameter than those tested during the present study. ‘The values for “n” show 

a slight decrease for pipes of large diameter — the over- -all change between — 

1 and 10-ft diameters is from 0.025 to 0. 023, whereas “f” values decline a 

sharply from 0.117 for 1-ft-diam pipe to 0. 046 for 10-ft- -diam pipe. Parallel 
if curves for the 25 per cent paved pipe are shown, in which “n” equals approxi- 

on ay 0.021 and “f” varies from 0.053 to 0.038 for pipes 4 and 8 ft in diame- 


ter, respectively. No « carves are ‘shown for 50 per cent paved pipe because 
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on 
aad ‘Sead It was recognized that several factors might have affected the accuracy of 
| i” the test results. Perhaps the most important of these factors involved the ee 


of calibrating larger orifice meters against several 


e condition throughout the course of the study. In addition, the random elgg 
ao) of the test pipes prevented an exact study of flow characteristics. For ex- _ 


a Z in. farther into the pipes than did the rest. Such inequalities within the test 
- had no effect | on or roughness coefficients, 
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a pa % ample, | a 6-ft-long straight-edge | laid along the inside of the respective pipes 


a smaller ones. : However, every known precaution was taken to ensure that all 
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but their presence the nonuniform of flow that 
sometimes was evident in the traverses. 


CONCLUSIONS 


coefficient Aah the range of experimental discharges 
reproduced for study. The following values of Manning’s coefficient “n” were 


_ obtained from tests on standard and — -invert corrugated metal pipes 
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REVISED COMPUTATIONOFA 
‘VELOCITY HEAD WEIGHTED VALUE 
Joe M. Lara,! M. ASCE and Kenneth B. M. 
q 


revised procedure is proposed on the of a 
weighted value of the velocity head. An example of the procedure is illustrated — 


by a hydraulic of the a m using the ‘slope- area 


“This paper a suggested | in the 
now used by hydrographers to determine the discharge of a stream by the ce 
islope-area method. A simpler way of computing the weighted value of the 


ivelocity head in channels —r overbank or subsectional flow is sho 


r sub = 


The usual method for determining the weighted value of velocity head is by i? 


» CON 


where Ky, Ke Kp, and the conveyances for each of the 


lof Areas Ap and An respectively. 


Note: Discussion open until February i, To extend the 


written request must be filed with the Executive Secretary, ASCE. Paper 2149 nae 
part of the copyrighted Journal of the Hydraulics Division, Proceedings of the Amer- 

Society of Civil Engineers, Vol. 85, No. HY 9, » September, 

2. Hydr. Engr, Commissioner’ Office, ‘Bureau of Reclamation, D 4 
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is the cross-sectional 


the roughness coefficient. 
The above method assumes that the velocity within each subsection 43,6 or 
_- Nis constant and that the discharge o of € each subsection is is Proportional to: its m 
computation of K. As As an example, , for K = 230,000, KS = 12,200 _— 000, 000, 000 
as.: the following procedure is proposed for co computing the weighted value o 
7 oe velocity head by a more convenient method. The proposed procedure is 4d y 
computations. (1) This procedure can be best illustrated by an example using 
A the slope-area data collected at Sections 2 and 3 of a stream. : A table of the — 


‘The equation for a@ reveals that astronomical values can result in the e) 
7 similar in principle to that used by ; the Corps of Engineers in backwater — 
‘sample computations and the - descriptions of the determination of the value in 


each column are § shown on the 


Section designation. Sections 2 and 3 are the 
upstream crc cross sections, 
Observed water surface elevations at each 


Subdivisional conveyances of each cross 


for example on line 1 


Conveyance or K is defined in terms of Manning’s formula as: | 
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fe — Column 7 divided by Column 3 q 


; tu Column 8 squared times the corresponding 

hted velocit head value com- 


pve puted by dividing the sum n of Column 9 by the 


| 
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Columns 15 and 16. When this value agrees 

a: reasonably with the assumed discharge (Column 

the are completed. If the re- 


‘To compare this , suggested procedure with the usual one, the following — 
shows h how a would be derived for 3 in the 
000)3 + (3,670)3 + 


, 490,000,000 


This value multiplied by results in 0. 466, the weighted 

‘velocity head for Section 3, the same answer as obtained in the example. — oe 

addition, another like computation must be made for Section2, 

oa The above suggested procedure provides additional information not — 

available from the usual method and that is—the velocities and discharges 

® occurring in each subdivision are readily computed. True, they may be — 

theoretical, , but they indicate how much of the flow is conveyed by « each seg- it =. 


ment. This enon is particularly useful for the Sae of back- 
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EDDY DIFFUSION IN HOMOGENEOUS TURBULENCE 


Taylor’: theory of diffusion by continuous movements and the Einstein 
equation of diffusion were applied to eddy diffusion of particles in a two a any 
dimensional field of homogeneous turbulence produced by broad shallow chan sg 
nel with extreme bottom roughness. A simple method of determining [ie 
Lagrangian: eddy size, eddy diffusivity, turbulence intensity, and mixing length 

was devised. "Application of this technique provided experimental verification _ 
of the Kolmogoroff similarity principle which relates eddy diffusion, the rate = 
of energy dissipation, and the Lagrangian eddy scale. _ ss of “a 


An interchange of fluid between two neighboring sla tn a turbulent field 
necessarily results in the simultaneous interchange of every characteristic — 
of the fluid masses involved. For example, if a fluid mass labeled witha a 


boundary into a zone which is free of contamination, it must, in order to sat- 7 
isfy continuity, be replaced by a fluid mass which contains no fluorescein. end 
The immediate result is to increase the average | concentration of t the dye on fe = 


one side of the boundary and reduce it on the other, each region acquiring the — 
properties of the other. The process, if allowed to continue indefinitely will, 


of course, result in producing a uniform concentration of dye rae both 


Ina turbulent. ‘field, the masses s involved in the interchange are called 
* “eddies? and the process is known as eddy diffusion because of f its si similarity 
to molecular diffusion. ‘The ‘primary difference between the two processes _ 


lies in analogy d. drawn between | discrete or molecules, moving 
| Note: Discussion open 1, 1960, To extend the closing date o one month, 
Written request must be filed with the Executive Secretary, ASCE. Paper 2150 is me 
E % _ part of the copyrighted Journal of the Hydraulics Division, Proceedings of the ihe ae 


_ American Society of Civil Engineers, Vol. 85, No. HY 9, September, 1959. 


. Asst. Prof. of Civ. ‘Eng., Univ. of Calif., Calif. 
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space, on the one and the rather defined lumps of fluid which 


> diffusion. Ina turbulent regime, the effect of molecular diffusion i 

ae: The rate of change in 1 concentration tie eddy diffusion is determined — 
_ fundamentally by two independent variables: the carrying capacity of np ; 
flow and the extent of difference between the concentrations on 
either | side of the arbitrary boundary. ' The latter variable can be expressed 

mathematically or evaluated experimentally as the lateral rate of change of : 

_ Concentration. The former is defined as the lateral coefficient of eddy dif- 


as discrete entities which behave in a manner analogous: 
to particles experiencing Brownian the Einstein(1) is 


where Dz is the lateral coefficient of eddy diffusion, t is time the oz is the # 

standard deviation of the lateral displacement from the mean position of a 4 


statistically significant number of eddies. ih 1 a turbulent stream, eddy mo- § 
_ tions may y be traced by introducing particles ofa foreign substance which can 
are _ be readily identified and the motions of which are — to those of the 
 Fick’s laws of diffusion, which were derived from with equations 
bony heat conduction, | have also been commonly applied to eddy diffusion. Fick@ 


and t ms Dz are as previously defined » has | been especially popular. wa a 
_ should be recognized, however, that when D, is not constant but is a function 
of t or z, solution of (2) is difficult and in many cases cannot be accomplishé 


es fen ‘The mixing length concept of fluid mixing, usually credited to Prandtl, .g 
ecognizes an analogy between fluid turbulence and the “ ‘mean free path” con- 


y 
ept of gas dynamics. This results sina 1 definition of diffusion — 


or ‘mixing | ength is the ro e root mean 

‘square velocity fluctuation in n the z- direction. a This relation for Dy possessé 


the ed of f incorporation of a property of the yes the velocity fluctua- 
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4 the ‘result that the mixing length idea way to statistical 
descriptions of mixing phenomena. 
‘The modern day statistical theories of turbulence stemmed largely fro 
Taylor’s theory of diffusion by continuous movements presented in 1921. - ha 
Taylor the of particles of fluid by means of the 


scheme and Re, , the velocity fluctuation correlation coefficient is defined ma 


§ failure of the Fickian laws of diffusion to account for the extremely re 


range in values of the diffusion coefficient observed in atmospheric turbulence. ie 
He introduced the concept of “neighbor distribution” which for an average 


separation of pairs, Lk results in the neighbor diffusion equation i 
ih 


where q is ‘the mean of unit length of L, the neighbor 
separation, and a * a diffusion coefficient which relates to the scale of -* 


in which e was taken. ran a constant. This relation, although presented with- 

| out basis in theory, ! has Ss been designated 2 as the ~“Four~- 45) a, Law.” It has 

’ been applied to oceanographic diffusion by y Richardson(4,5 by and others(6) v with ie. 

| varying degrees of success. The exponent four thirds was mesrutsay by some in- 

vestigators as a universal number ‘despite the > lack of of Satisfactory experimen- 

tal support and the absence of a theoretical basis for this conclusion. or. 

_ The four-thirds relationship remained without prospect of theoretical sup- 
port until 1941 when Kolmogoroff(”, 8) introduced a somewhat new theory of 
ae _ Gj the spectrum of turbulence. Kolmogoroff’s s , development is based on the fund. 
t mean § mental concept that eddying motions are characterized by a wide range of 
possesse length scales. Energy « of f eddy motion is then conceived tc to be passed down th 
fluctua- spectrum through a succession of eddies of smaller size » ultimately being dis Say 

sipated by viscous action as heat. In this process, as the scale is pope ae a 
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y fluctuation and é is the time displace- 
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les, Richardson suggested that 
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‘The 


types of turbulence. By such eddies are isotropic. = 
“a Kolmogoroff’s analysis is built on two hypotheses, both applicable to mo- 


mt The small-scale components of the motion depend only on the viscosity, 
and t mean rate of unit mass of 


to the scales oft the ‘motion, thus 


can be replaced with D. ‘By mare results 


literature fails to reveal an confirmation of the 


tal Tech j 


= 


: ipation per unit mass of fluid, | and the coefficient of eddy ate | 


ns homogeneous turbulence in three dimensions is | virtually im- 


large conduit probably come clowest to the ideal requirements, but ‘boundary 


effects severely restrict the scales of turbulence which may be examined. 


_ Wind or water tunnel grids may be expected to produce homogensity of turbu- 
lence in the field stream but only in planes parallel to the grids. ‘Because the 
- scale of turbulence diminishes with distance from the grid, this form of turbu 
_ lence generator is not suited to a study which attempts to relate eddy scales 


and eddy diffusion at a constant rate of energy —<— 
‘If the requirement of three dimensional homogeneity is relaxed and a two 


dimensional field is substituted, the experimental techniques are considerably 


‘simplified. Al broad, ‘open channel 0 of suf sufficient length and bottom 
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roughness to minimize end effects meets the requir ft 
homogeneity very satisfactorily provided the diffusion phatalens | i 
strained to a single plane, i.e. the surface of the stream. 
_§} The experimental channel constructed for the studies was 4 feet wide, a 
scosity, inches deep and 38 feet in length (Figure 1) and was supported on leveling 
pea screws which permitted adjustment to uniform slopes ranging from 0 to 0.01. 
7 Two one-quarter horsepower propellor pumps were capable of delivering 
depend § flows sup | to 1000 gpm. to the e stilling basin at the head of the flume. — _ A by-pass aie 
id. a valve between the stilling basin and the return flume facilitate edischarge ‘ 
nly be om Bottom roughness for the experimental channel consisted of expanded — Core 
aa “metal meshes as shown in Figure 2 which were stapled to the plywood bottom > 
eS tel of the channel. Three mesh sizes designated by the nominal dimensions — - 
the minor diagonal of the diamond pattern as" 5, 1.0 and 1.5 were 
4 _ The requirement that the diffusion phenomenon be observed only ina plane — 
= where two-dimensional homogeneity could be achieved dictated the use of = 4 
§ floating particles, the motions of which were influenced by turbulence near - 
© the stream surface. Use of discrete particle motions to statistically charac-_ * 
@ terize the turbulent regime is in accord with the theories of Taylor(2) and an 
® Richardson. (3) Also, a a simple direct method for statistical interpretation of 
particle motion in terms of eddy diffusion is available through Einstein’s 
ionships The particles employed in the studies were ‘polyethylene discs approxi-— 
lya mately one- eighth of an inch in diameter and one-sixteenth of an inch thick. _ 
“ahs Polyethylene has : a specific gravity of approximately | 0.975 so that, when fully s 
wet, discs of this material just barely float. The buoyancy and size of the . 
discs used in these investigations were ssbecathe satisfactory to prevent 
their being entrained in the flow and to insure that their motion was repre- 2 ve" 
sentative of the upper er few millimeters of the flowing stream. Particles 
; released singly through a small glass funnel suspended about one- eighth of an oe 
po inch ¢ over the water surface from a point gage assembly w which could be 45 a £ 


owever, 
ry under 


_ In accordance with diffusion theory, the position of each particle BAN “a 
fixed ata certain time following the instant 0 of release. Each point is pal 
projected toa single a axis. ‘The distribution of particle positions along oll 
axis can be defined then in terms of the standard deviation from the mean 
. | position. The procedure used in the investigations reported herein differed — 3 
§ slightly in that particles were all intercepted from the flow at a fixed dis- = 
tance downstream from the point of release. However, since the theoretical - Pet 
pattern would be symmetrical about this | arbitrarily selected axis, the total Y = ue 
joundary experiences, i. i.e., the total accumulated path length of all particles, would be ~ ae 
Gidentical by either technique. It is, therefore, reasonable to conclude that the 
nined. 
; of tu bee standard deviation calculated from an array of particles along a fixed line of a a 
interception a distance, %, from the source is substantially equal to that of a iS. 
pattern whose centroid is also x distant from the source. = ~~ a 
Particles were collected from the surface of the experimental channel by aor 


cause the 
us 


§ inserting a scoop of ordinary galvanized window screen into the top one- 
nd ~~ quarter inch of the flow. In order to simplify the statistical analysis of the a ag 


nside 
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etal sheets spaced one centimeter apart across the full width of the channel. ie 4 . 
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a _ In order to reduce the time required for determination of the standard — : 
eviation of particle position, ordinarily a very time consuming operation, 

imple graphical procedure was adopted. It consisted of accumulating the 
- particles captured in the collector from one side of the dispersion ome 
a expressing them as per cent of total numbers, and plotting on arithmetic 
_ probability paper at intervals corresponding to the spacing of cells. The re- 

; sulting plot was invariably of the typical Gaussian form; that is, it could be 
closely fitted with a straight line as shown in Figure 3. ” The standard devia- 
_ tion was then estimated graphically from the slope of the line. A comparison 
of this technique with the formal computation technique for determination of a 


a _ indicated the standard deviation in the error of the graphical method was 


OTAL 


T 


a TICLES ACCUMULATED - PERCENT 


- clusive or incomplete results and thus producing maximum benefit from the 


Prior to the actual investigation, the performance of the experimental 
an ana was verified in order to Goterarine whether or not two dimensional 


homogeneity was actually attained. Determinations « of at fixed distance 
from the source over the width and along the length of the channel indicated 

a ue that except for points very close to the channel walls, say less than : inches 

from the walls and within 5 or 6 feet of the entrance or exit channel, the tur- 


from the of release thus experimentally, the relationship be- 

‘ a Og and x. According to the basic Einstein diffusion equation (Equation 

oe om is $ directly proportional to x1/ 2 when the coefficient of diffusion and 
stream velocity are constant. Both the theoretical analyses offered by ; os 
Taylor( 2) and by Richardson(3) agree | with this relation for large values ae 


prensa experimentally would be wedge-shaped near r the source and would 
_ gradually approach a parabolic pattern as the distance from the source is in- 
— 4). Also, the diffusion coefficient, which may be considered 


> 
— 
d 
j ing so great that the graphical technique was fully justified. Moreover, the 
method facilitated interpretation of experimental results while the experiment 
— 
— 
7 q sy observations were confined to a smaller area about 2 feet wide and 15 feet 4 
| in the central portion of the — | 
Using the methods described above, it was possible to determine the 
— 
— 
— 
5 @ _ to xX, ultimately attaining a constant and maximum value distant Irom the = 
If the Einstein equation is modified by substitution of dx/u for dt, the 
diffusion coefficient becomes related to distance by 
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ao Dz 2 may be expected to be constant t and, in particular, how 0, varies 
= Such a as those shown 
2% col 
a) in in which the er n represents the slope of the curve and m is a 
9 "Differentiating | (1 1) with respect to x produces the general equation mm for 
the diffusion coefficient at any distance x, fromthe source 


DISTANCE 


in which the and n characterize the tangent to the o 2 -x curve 


be regarded as a crude measure of the scale of the 


ges 
fusion then by replacing x in equation: 12 by its equivalent in 
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SMOOTH CHANNEL BOTTOM 


5 0, = sim 

if n=l, the special case of constant diffusion results. 2, 

and o, are directly proportional to x, the distance from the source, 
-_ nee as suggested by Taylor and Richardson for small values of x. And finally, if 
i, es cn = 3, there results the special case of the diffusion coefficient proportional ae, ae 
to the four -thirds | power of the diffusion scale (co, e.astatement ofthe 
Z/? 
ty "similarity principle of turbulent for a rate of 
oe In Figure 4 the value of n is seen to range iba 1.0 to 1.59 ik ae 
x and on the Reynolds number, . . There were indications in these 
and all similar series that at sufficiently large values of x the value of 


‘It was apparent from prelim 
fication of a “four- -thirds” law or the ‘similarity principle would not be eae 
ble by direct observation of the dispersion pattern in laboratory scale chan- 
‘The si ib-range of eddy sizes to which | th imilarity would apply 


tus indicates by its slope the regio 
te 
a 
4 
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stream through a a fine glass the shape of dispersion pattern 
_ closer to the source could be examined. Such dye streams were photographed | 
with time exposures on high speed film under ultra violet Me icic A baths q 
- Assuming that the dispersion pattern produced on the _ photographs defined 4 
limits: of + 20z z, it was noted, as shown in Figure 5, that the value of n 2f 
a4 frequently exceeded 2.0 and occasionally approximated the value of 3.0 
theoretically required for conclusive verification of Kolmogoroff’s relation. — 
_ It is clear that direct corroboration in the laboratory of the applicability of | — 
: equation 9 by inspection of individual dispersion patterns will require much | 
finer measurements than were possible in these studies. 
‘ However, the relationship between Dz, the eddy diffusion coef ficient, a and 
q ‘ La, the mean eddy size, can be "gutirunase through a a ‘comparison of the 3 
i oe properties of several different turbulence fields. In order to accomplish | this 
; A comparison it was necessary to relate characteristics of the dispersion pat- 
tern, Oz, tc to the size of which brought about dispersion. 
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In order hand side of (14) it is essential to know 
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uation of the form 
Re 


Substituting into (14), and to with hx/i, produces med 
Payot. 


a 


With the aid of equation on (16) and an experimental knowledge of Ce, % and u, 


vy ‘derivable from the basic equation. it will be noted that the ‘shape of the 
-x curve is determined ed only by the eddy scale reciprocal, the velo- 
ba city fluctuation being only a a translational term which positions the curve gt 
- along the vertical axis. Thus, k/i can be determined without knowing wl2 0% 


by ‘Simply ‘fitting a a curve of relative 2 versus x x to experimental data over 
ia suitable range of x. x. In fact, one nnet only know two pairs of values of of i 
and x to obtain an estimate of eddy scale. Once k/ii has been established, 
tes w'2 can be e calculated from Equation | 16. - Since the velocity is ordinarily 3 7 
known, can be permitting of the eddy 


have hee h sma er even: han size of par ies oved 
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Raichlen(17) show that the Bé curve may be closely approximated by an 
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vestigation, a | group | of normalized o 2 -x cu were prepared, samples of 


- which are shown in Figure | 6. The ordinates for each of thet: curves are the 


values « of o2 relative to the value at x = ‘The abscissas. cover the normal 
"experimental range from 1 to 20 feet. Each curve is labeled with the para- 


meter k/i. After plotting o2-x data for a given run, an overlay carrying , 3 
_ the normalized curves at the same scale was superimposed and the appropri-_ 


_ ate value of k/i was selected. ‘Figure 7 shows a typical set of experimental q 


= fitted with the general curve. As might be expected, the formulated ie 


A curve agrees with the experimental data fairly well within the range of 1 bail 


feet from the source in the region close to the 
While it was not the relationship between eddy size 
and eddy diffusion for an individual dispersion pattern it appeared feasible to 
attempt relating the ultimate coefficient D, and the Lagrangian average 
eddy size, La obtained from a number of independent experiments. Although 


_ the average eddy size, characterized by the area under the Ré - - € curve, is 
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the so-called “Four - ‘Thirds Law.” Also, it is of some significance to ob- 


= U is the mean stream velocity, 

is the slope of the energy grade line. 
| es In the experimental channel it was virtually impossible to attain precise a ioe 
f control over either 1 E or Ly since br both | quantities were affected by changes" 
¢ in » slope, velocity, and channel roughness. To examine the independent vari- 
ability of of E with Dz (0) for example, i it was necessary to relate Dz’ (0) and 
La empirically fora number of of channel : slopes and roughnesses, cross-plot 

“ther resulting curves | for constant Ly! and compute the associated values of E. 


Figure 8 shows typical examples of Dz (co) - ~ . La curves. of the type » usually Ae 
- associated with the Richardson relation (Equation 7) in which E is not held — e 
constant. _ ‘It is interesting to note that each set of plotted points can be fairly Sih 
7 well represented with a line of four-thirds slope, an apparent corroboration i. 


_ serve that the vertical position of the curves is affected significantly by the - 


. 3 Boren of the channel bottom but apparently is little influenced by the size of 


bottom roughness. - ‘The fact that Dz (ce) : is dependent o on slope, i. e. is propor - Oe 


“tional to the rate of energy dissipation, generally with the similarity 


concept. 
7 Wve. computation of E for constant values of La also required empirically Nae 
relating the mean channel velocity, to be and in turn selecting Se 

a channel calibration curve of Se versus channel slope, and channel velocity. 


4 


of E determined from equation 22- and corresponding t to eddy scales og 


4 


of 0. 2, 0.5, and 1.0 feet were plotted against the appropriate values of Dz (co) ae , 
i? Despite the fact that each plotted point in Figure 9 is supported by 50 to ae 
100 separate determinations of the standard deviation in the dispersion pat- 


tern there is considerable scatter in the plots. This must be attributed, 


ne plots : appears to ag agree a well with that expected from theory. ‘taf Rk 
_ Figure 10 shows the coefficient of lateral edd A diffusion as a function of ie 
Lagrangian eddy size at constant E (0.05 ft2/secs). _ The trend line in the rus 
figure is seen to have a shape similar to those of Figure 8 approximating - S 
four-thirds slope at large eddy sizes and tending to flatten at eae .. _ 
Again, the experimental data appear to support tl the theory. 
A further test of the validity of the similarity principle is by 
Toon values of E1/3 L,4/3 against D, without regard to the slope of the a By 

_ channel bottom or the condition of bottom roughness. ‘Figure 11 ‘summarizes — 
the data representing 105 separate experiments in which-U, and bottom 
Toughness are all independent variables. It is especially significant to note — ay 
- the plotted points can be fitted with a straight line on a unit slope. T = ae 


the line which best fits the data of these experiments 


i where re Dz, is the coefficient of ‘the lateral 1 eddy diffusion, E is is the e rate Be 
energy ‘dissipation per unit ‘mass, and La is the Lagrangian eddy s size. 
constant is taken as representative of all degrees of bottom roughness pro- eS Ah: 
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The. fact that all the data the series ant to pre prepare llare 
_ rather equally distributed about the best fit line indicates that the scatter is ze 
_ probably due to experimental errors alone. It may also be concluded that a as 
_ independent variables are probably accounted for in the quantity E 1/ 314/38 | 


ort that the experimental range of the variable was not sufficient to cause i. ie) 
separation in the plotted points. The latter conclusion is probably applicable 
4 in the case of bottom roughness, the absolute value of which was not greatly 
- affected by a a change in mesh size. There is some indication from | exploratory at 
BS tests w: with a | smooth channel bottom that the vertical position of the line of best ‘ 


An diffusion | of stream 3 “4 


‘we definition of the coefficient of eddy diffusion and the development of a tech- 

— nique Satisfactory for its determination. There are many ways in which Dz 

e can be defined but most of these require the separate evaluation o of quantities i 

b such as the turbulence intensity, the mixing length, or the eddy scale in cr 
4 - that actual values of the eddy diffusion coefficient it may be calculated, Be- 

cause of difficulties involved in the physical measurement of these quantities: 


_ there has been a ne for investigators to resort to empirical or : * 
4 


5 

oF 


|  0900327,05" 
0.00163 


(0.00055,0.5" 
0.00055, 


fo Sey D. = 0.0136 


= 


NT OF LATERAL EDDY DIFFUSION , 0 


FICIE 


| 
EC 
‘AS A FUNCTION OF SCALE AND 
OF ENERGY DISSIPATION 


tt 
3 
— 
— 
| 
— 
2.28 
TAS 


and constants of proportionality and appro- 
a priate exponents. Often the empirical equation is only a convenient device 


intuition and astute observation, hit upon a caebiacian which comes close to — 4 
— the truth but which is not quite complete. ‘This appears to have been the case 
with the semi-empirical relationship first proposed by Richardson and which 


has been referred to as the * Four-Thirds Law. Rees 2 


‘ ree Richardson observed that the coefficient of eddy diffusion was nearly pro - | 


‘portional to the four-thirds power of some measure of the scale of the phe- | 
fe a nomenon causing diffusion. He proposed a theory of “neighbor diffusion” _ a 
ie which involved the statistical averaging of “neighbor separations” of adjacent | 
‘particles. _ The mean neighbor separation he took to be a measure of “scale.” 
«Unt jrtunately, Richardson’s neighbor aiffusion theory is difficult to apply, and 
" except for a few studies on ocean surfaces, it has been little used, oa a 


_ Pearson(6,1 18) has summarized the relationship ‘between the coefficient of ; 


Br bray of many independent investigators, including h himself, fo 

variety of scales. Figure 12isa typical plot of the type Pearson p 
including all of the data reported by him together with some observations ~ 
woe collected during the investigations which are reported herein. It is possible 

to see from the figure why one might be inclined to regard Richardson’s re- 

c ation as a universal law. A line of four-thirds slope fairly well establishes" 
the over- -all trend of data | from scales of about 0.1 foot to scales of more than 
: re 1000 feet. . The considerable scatter in the separate observations can probably 

be attributed, in part, to the variety of definitions used for “scale.” These 

definitions included mean neighbor separation of particles, one- quarter of the 
width of the dispersion pattern (approximately = — oz), the Lagrangian | eddy 
ize, and some arbitrary dimensions such as width or depth of channel. oe . 
Fs It has been shown theoretically according to the similarity concept of e 
q Kolmogoroff that the fundamental relation should include a third variable, 
s  _—séoAnthe: rate of energy dissipation per unit mass. It may be argued that much of | 
Ls: a the spread in the observations of Figure 12 could have been eliminated had 

ty PS the rate of energy dissipation been considered in each case. The most ap- 

__—s propriate abscissa for a plot of a supposed universal law would be | EV 3 L4/% 3 
with a consistent definition for D, E, and L as in Figure 11. ae eae -* 

a It is also apparent from consideration of the shape of the dispersion- sit 


enon, but that it is limited by the maximum size of eddy motions which may 
exist in the turbulent regime. The size of these limiting eddies i: is not a 


accordance with the similarity principle, there 1 will be a sub- ~range of eddy 4 : 


sizes to which the general diffusion law will apply. _ Within this sub- ~range it~ 
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_‘The significance of the of a ‘diffusion coefficient can best 7 
be seen by considering the problem of waste assimilation by a receiving water, 
gaya a large coastal bay. If the waste is discharged ata a point on the surface 
the pattern of dispersion may at first appear to increase rapidly, asin 
Figure 13, fanning outward as the waste diffuses laterally while being carri 


= by the current. In the region near the source the similarity sores 


4 will no longer be cau tee may y continue to be influenced by eddy size i 
_ but to a lesser degree, and the dispersion curve will gradgally change guage, 3 


i will have become involved in the retrain pattern and the eddy « diffusion 


ie coefficient will cease to increase. Beyond this point: the ‘dispersion curve will B 


parabolic in in n shape, a constant 


‘The Kolmogoroff similarity concept leads" to the eoncleaion that 
a the coefficient of eddy diffusion in homogeneous turbulent flow is prepeetioniil 
Pe. to the one-third power of the rate of energy dissipation per unit mass and the 
four-thirds power of the scale of eddies participating in diffusion. Using a : 
open channel with a variety of bottom roughness to produce a two- = 
dimensional homogeneous turbulence, the e experimental validity of the concept 
examined. The following conclusions resulted from the investigation: j 


1. In laboratory turbulence the ‘Sub-range of eddies to which the concept 


The upper. size _ 
_ limit in the open channel studies reported was probably less than one centi- 
_ Examination of particle and dye dispersion nas we close to ‘the 
source indicated the of the sub- “range. 


mit of the critical sub-range conform closely to the peed a of Taylor all 
Richardson. Close to the source the relationship between the standard devia- 
tion of the dispersion pattern and the distance from the source is almost 4, 
linear. ‘Far from the ‘source the dispersion pattern becomes parabolic it in 
; shape. In the transition range the pattern can be described statistically by 
assuming an form for Taylor’s s Lagrangian velocity correlation 

2. ‘For each spectrum there is is a maximum value for the eddy 

diffusion coefficient which is associated with the limiting eddy size. The 
Timiting eddy size is directly proportional to the statistical mean eddy size. 

_ 4, The Kolmogoroff similarity principle adequately describes the relation- 
ship between the ultimate lateral eddy diffusion coefficient t Dy the 


Lagrangian mean eddy size la, and the rate of energy d ‘dissipation per unit [oe 


mass E, in statistically homogeneous s flow i in broad a _ mera ie 4 
equation 
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WARREN B. eeaianes M. ASCE.. —The writer is indebted ae 
Fred W. Blaisdell and Harold W. Humphreys for their interest in preparing a. 
discussion and i in offering fo for consideration the drop inlet type of structure. ee 
- - Since much of the discussion centers around the definition of a siphon, ot 
might be well to analyze the two definitions given. The first, to which _ 
writer subscribes, ‘defines a siphon “ “as a spillway which requires | that the “a 
“- “liquid be lifted above the headpool level.” This definition is supported a 
Russell,(1) Addison,(2) and Rouse.(3) The Bureau of Reclamation so classi- 
| fies such a structure, though this position admittedly is based on i con- 
a _ The second Sediaition, which the euthors suggest as being a appropriate, is 
‘ fa closed conduit, a portion of which lies above the hydraulic grade line.” 
| Three Bureau ‘structures are named as qualifying as siphons, the Heart Butte 
| and Shadehill morning- glory spillways in the Dakotas, and the Hungry Horse 
morning - glory in Montana. The writer patieves none of these operate as et 
siphons, even by the authors ’ definition. All three | spillways were studied ex ¥ pt 
_| tensively by model experiments. In general it can be said that certain por- — ' 
| tions: of the profiles showed negative pressures for ‘the unsubmerged flow — 


received radial or near- -radial flow, and not between certain elevations for oo 
full circumference. Thus, the low pressures were similar to those on an ogee — ; 
shape, and were related to tendencies for flow separation from the boundaries _ 
because of imperfect profiles. Since the upper nappe surfaces were essen-_ 


tially at atmospheric -dpearpeg no siphon head could exist because the flow a “nt 
| wasoftheopenchamneltype; 


4 ‘The possibilities of extremely low pressures occurring with 1 submerged 


| flow conditions were indeed recognized. At Heart Butte a combination of nine __ ig 
| air inlets and deflectors was installed 11.5 feet below the crest and an aerated aA 
flow deflector placed on the crown side of the shaft 31.5 feet below the crest. | F 
Further, the tunnel diameter diverged from 11 to 14 feet in diameter om ong “— 
the bend. Pressure tests at submerged flow showed negative pressures no - 
greater than about -5 feet of water in certain areas, with no indication that : 
low pressures existed around the entire circumference, as would have been vn 
the case if a closed conduit siphon head had been acting. 
‘@ There was less opportunity for low pressures in the Shadehill | spillway be- 
cause of the large spillwa way and tunnel diameters and the short drop of 33 feet 
a. Proc. Paper 1807, October, 1958, by Warren B. ‘McBirney. 
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exceeding - -3 feet of water 
- occurred for weir-flow discharges less than 60 percent of design capacity, 
and all measured pressures were well above atmospheric at capacity flow. _ 4 _ 
_ Even so, a flow deflector 1 was installed on the crown side of the shaft 14 feet gy “let 
below the crest to insure open channel flow from the deflector downstream to . 


ae The Hungry Horse morning-glory was studied in one of the most thorough 
investigations ever conducted by the Bureau. The negative pressures on the e 
ae profile were deliberately created to increase the discharge coefficient and a 
ifn were restricted to areas below portions of the crest receiving radial or near- 
. flow. Other negative pressures were relieved by an air vent. An ar a 
_ intake structure of 36 square feet cross- -section was installed on the crown a 
of the shaft 31 feet below the crest. conduit had sufficient air 
capacity to prevent the shaft and tunnel from flowing full for any flow or crest 
_ thought possible at the site. The writer as well as the research engi- 
ae _iheers responsible for these tests believe there is no possibility of the all 
Horse spillway and tunnel ever flowing full. SA aj 
None of these structures in model or prototype has shown any evidence of 
_ siphonic action in the sense that a full conduit was acting to influence the rate | 
if 3 of flow at higher elevations in the shaft or tunnel. No priming or depriming _ 
_: ia aids were built since the structures were not envisioned as siphons and every | 
effort to prevent siphon action was made. ‘The air demand of | each structure | 
was more than satisfied with oversized conduits to prevent unsteady rates of : 
; ae . This is not the type of partialization desired by the writer. ee 
_ With the above description of the operation of these spillways completed, — 
_ the writer is now prepared to admit that a morning-glory or drop inlet spill- = 
_way flowing submerged with an unaerated shaft of sufficient drop in elevation 3 
will exhibit flow characteristics akin to those in the ‘downstream leg of a sad- 
_ dle or volute siphon. It is doubtful, however, whether the engineer who pro- | 
duces such a cavitation structure would derive much credit for his design. a 
a _ The head-discharge curve in Figure A of the discussion, for the drop | inlet 
closed conduit spillway, is very similar to | a typical morning- glory head- — — 
si discharge curve. Placing an antivortex flat plate over an inlet 3 inches in 
Ke diameter may appear to offer control of the vortex cand a uniform proportion- 3 
ing of air and water. Trial of a similar scheme in a model of the Trinity Dam ow 
morning - glory, northern California, model size 24 inches in diameter at the € 3. "4 
crest, brought near ~disastrous results to the clear plastic of tl the shaft. Air : 
was sucked under the plate intermittently instead of steadily, pressures in the 
shaft varied rapidly between -10 and -30 feet of water (prototype), the forebay 
: supplying the structure ¥ was soon § set in rhythmic motion, and the the ‘model vi- 4 
brated violently. action in the prototype would cause serious damage 


_ ‘The: writer suggests that the authors employ a high frequency response bong 


pressure cell and t to a piezometer centered in an antivortex 


Near th the end of their discussion the authors ewes a linear + multiplication 
factor to compare operation of of the Bureau saddle siphons with the drop inlet © 
spillways. Such a by this means 1s not believed to be valid for 
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both types of structures to be siphons, scale relation- 


ships for siphons are still very much in doubt and are probably not linear. a “i 


S r 2. The authors compare by the multiplication factor a head on the drop =a 
_ let which “caused the Illinois spillway to prime ” to the heads for maximum > 
— flows in the standard and proposed eg ~The latter heads were 


‘any information on the effects of changing the rates of rise in the forebays of 
Tespective spillways. ‘The Bureau test setup was not flexible enough to =. 


been able to do so. , Findings are that the lower the rate of rise for a given 


1 discharge the lower will be the head required to couse the : siphon | to prime. 


_ at lower heads for lower ratesofrise. = 


The authors state that extending skirts from the antivortex 
tT below the weir crest modifies the structure to the form of a vortex siphon — 3 


spillway and that the head-discharge curve then obtained is similar to those 


tures and that only one of them should be classified as a 2 aiken Gat 


_ The writer does not doubt that the drop inlet closed conduit spillway has 


x uses , but in its current form it does not appear to fit oe = 
the Bureau of Reclamation for an emergency ee: : 


. “Hydraulics,” George E. Russell, p 213, Henry Holt and Company, 1942. n 


A Treatise on Applied Hydraulics,” j 
Wiley and Sons, Inc., 1945. me, hee 


‘Blementary Mechanics of Fluids,” 


— 
_ started to prime with about 0.13 foot of head and the proposed siphon re- 
‘3 | quired about 0.015 foot of head. Both heads are less than the 0.17 foot shown — &g ae HE 
the authors for the Illinois spillway. 
& — 
— 
— 
rate 
| 
ery — 
— 
veal — 
sad- | 
inlet | i 
erbert Addison, pp 200-202, John ss 
the “wi Wiley and 
— 
| 


eri 


a 
— 
— 


Rita A 

OF LAKE ONTARIO® 


MAX A. ASCE. 
_ methods of computing -evaporation— mass transfer, water budget and alll : 
4 budget — -only the first two of which he has applied to Lake Ontario. . In addition 
to the methods cited by the author, there is the evaporation-pan approach as 
i poe as several others which constitute combinations of the four. It should > i 


equations which have been derived by numerous investigators. 


_ Considering, first, the mass-transfer approach, it will be noted that the - 
| author's equation (Eq. 7) relates evaporation to ) the vertical gradient of hori- 
zontal wind speed over the lake and that the empirical constant was derived 7 
- from data collected at Lake Hefner, Okla. (1) Although one might expect to ob re 
tain satisfactory | estimates of normal monthly wind speed over the water from 
— anemometers, corresponding estimates of the vertical gradients are a 1 q 


_ likely to be appreciably in error. In the Weather Bureau’s analysis of the 
_ Lake Hefner ' data, (2) several mass-transfer equations were derived empiri- a ; 


cally which require wi wind movement instead of wind gradient. One of a 


where E is the lake evaporation in inches per day, Uy, is the -meter wind = 4 
_ movement in miles per day, eg is the saturation vapor pressure at the tem- __ 
a of the water plone, Yard of mercury), and e eg is the 2- -meter 
vapor pressure over the water, also inches of mercury. The constant of this _ 
_ equation is based on records for 115 days when the water- r-budget « data were | a! 
_ judged to be accurate within 5 percent. The ‘correlation somal is 0.956— 
higher than the author cites for his equation. © 
__ Equation (A) has been applied to Lake Ontario using the ——re s data 
_ (Table 7) with the results shown in Table A. It will be noted that a (A) _ act 


sh 
_ 1. Mean annual free-water evaporation 4 


= Proc. Paper by Ira A. ‘Hunt, Jr. 
‘Chf. Research U. S. Weather Bureau, D. 
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| ~-Percent of annual class A pan evaporation occuring during the summer 


tember, 1959 


| 
= .150 


At Lake Hefner, 4-meter wind averaged 0.90 that at 8- meters. | get 


. + factor was applied to the data in the author' 8 ‘Table 7 NEY converted 


~de ‘b. The author's vapor-pressure difference ‘values are incorrectly aca 


‘ees as "millibars"., Strictly speaking, Eq. (A) requires the 2-meter nae 
: _ Vapor pressure over the water rather than the average for several 


The » conmpatatansidll we were based on. on records for the | period 1946- 55. In the 


_ course of deriving the free- water map, evaporation was computed from poe, 
radiation (or percent sun unshine), dewpoint, air temperature and wind for a. 


meteorological stations adjacent to Lake Ontario. The computed values 


would appear from these results that the 10-year mean from 
ak the lake is near two and one- -quarter feet— a value in close agreement with 
our ‘mass-transfer and the author’s water -budget results. bid 
‘Having essentially substantiated the author’s water- -budget value (2. 
by two independent approaches, it is indeed unfortunate that a question 


a must be raised concerning the value used for lake precipitation. Early 1 fewults tS 


from the island network on Lake Michigan have led the author to conclude that; Rc 
‘Lake Michigan and final isiieesel must await further analysis. Neverthe- 5 f 


less, indications ar are in accord with the author’ 


recor er, and annual precipitation 
P the island stations and comparative data for a group of lake-shore stations. 


4 There can be little - doubt that the islands (and the lake) receive less st nae 
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piscussion 
‘precipitation that would be indicated by an average of the shore stations. But 


i is this deficiency compensated by an excess during the winter months? The 
records from South Manitou (and possibly Beaver Light Station) indicate a. 
he (a _ The author concludes that at least part of the differences among the — P a 
data result from the fact that the larger islands produce an increase of oe 
i, cipitation. It would seem that any heating effect should be greatest at the 
: downwind side of a small island (and likely even beyond) whereas St. James _ 
‘ob , -and the Light Station are on the northern and southern portions of Beaver 
_ Island, respectively. Some of those concerned with the project have felt for a 


- some time that virtually all the differences could be explained by relative = 


‘ 
a In July 1958, anemometers were installed at four of the island ei a i 
Ke (gage-orifice height) which were read in the latter part of the following ae 3 
Mean wind speeds (mph) for the period were as 


James, Beaver Island (4.04 


Beaver Island Light Station 05 
Tle Aux Galets 02 


‘Table B. LAKE MICHIGAN PRECIPITATION, 
(SIX-YEAR MEANS, 1952-1958), ININCHES 
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den ‘St. James, Beaver island 


10.79 


years s (1956-58) for wi which the shore-station precipitation averaged 
ps 14.96, 13.22 and 28.18 for the summer, winter and enmel eget 
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Assuming that the average annual precipitation observed at a network of — 
eight stations along the shores of Lake Ontario is equivalent to that over r the 
a the value in Column 6 of the author’s Table 1 becomes 2.76 ft. The _— 
= evaporation then becomes 2.67 ft. instead of 2.21 ft. Since this 
revised value is appreciably higher than that derived by any of the other 
ethods, the writer has attempted to check the roe for a different period 


Precipitation (average of 8s ‘shore 2. 84 | 


Runoff (104% of 1934-53 as indicated by average 


Welland Canal) 


eee Outflow (St. Lawre ence River at Ogdensburg, N. Y. -) 38. 63 a 
Computed evaporation (1946- 55) 


Py ft is entirely possible that the / evaporation we was less in 1946-55 than for the 


coe longer period 1934-53, but ‘much of the difference (2.05 and 2.67 ft.) is believei 


thi 


to be erroneous. " Without belaboring possible sources of errors in water- f 


computations, it should be pointed out that bank- ~storage effects would | 


«ft. per year during the period 1946-55 and rose 0.11 ft. per year during 1934- 


oh ne «+58. _ Even 80, any realistic assessment of the reliability of the water -budget 


-— erapocation is not unlikely. A bias of only one percent in the inflow | _ out- 


items amounts to about 0.7 ft. over the lake. 
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‘Lakes, U. S. Weather Bureau Research Paper 38, May 1955. 


Wilson, W. T., Discussion of paper, Precipitation at Barrow. , Alaska, 
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han Re Trans. Am. Geophys. Union, Vol. 35, Pp- 206- 207, 


a ees In the light of Wilson’s(5) studies of wind vs. gage catch, it would now appear “3 a 
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_ PROBLEMS C CONCERNING USE OF LOW HEAD RADIAL GATES* 
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R. BOWMAN,* 1 ASCE.—Mr. Rhone’s paper r constitutes a significant 
contribution toward a better understanding of the hydraulic characteristics _ 
_ of tainter (radial) gates. 3 The author and the USBR should be congratulated + og 
_ for sharing with the profession the knowledge g gained through research | sup- — 
4 ported by public funds; in many instances, the results of basic research car- 
ried out by public agencies | has been labeled or 


Disc 


under a tainter gate, without recourse to a . model test program specifically e 
organized for this purpose. This is often the case for smaller | regulating a 
works whose total costs do not warrant expenditures for model tests. In 
_ few instances, it has been possible to check or modify analytical predictions ieee é 

by means of field measurements; more often than not, however, comprehen- 7 = 
sive field tests will not be feasible, with the result that the responsibility for 

‘discharge predictions will rest t solely with the designer. 


~~ * encountered as long as tainter gates have been in use. Where the a a 
sill is located at the weir crest, the rating curves for the gate will be mete 
dar to those of a vertical- lift gate, e, provided that the centerline of trunnions - 4 
is not appreciably higher than the centerline of the gate. Where the trunnions = a 
a _ are located at the elevation of headwater, or higher, the approach angle 6 a 

considerable influence on the value of the discharge coefficient, as _ 


shown in Figure 3. In cases in which the gate sill is located downstream ex a 


| ~ from the weir crest, a shifting control is encountered; The control point * 


a mm moves gradually : from the | gate sill, at small op openings, to the weir crest as 
full gate opening is approached. Cases in which the gate sill is located up- /. 

"stream from the crest are not often encountered, inasmuch as such a gate 
- setting ‘usually results in increased costs of associated structures which oe ‘ 
"more than offset any saving in gate costs. 
ee _ In the development of rating curves for a tainter gate whose sill is located . : 
_ at the weir crest, the writer, as well as many | of his predecessors, has tS 
_ utilized a method whose validity is largely confirmed by the data presented — 
a :. in Figure 4. _ When the approach angle @ lies generally between 75 and 105 ae 

degrees, the free. ‘discharge coefficient can be applied to all gate openings, 

= Proc. Paper 1 1935, February, 1959, by Thomas J. Rhone. os Be 
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with an average error of than percent. writer is acquai nted with a 


few such cases in which field measurements and pooticted: values differed by | shov 
a _ Although no simple and accurate method has yet been devised for comput- | 7 of di 
illus 

the weir crest, Mr. Rhone has presented information which may well serve 
as the basis for the effective utilization of a combination of methods. __ ll a 
e On the basis of analyses of a number of rating curves based on model ape | writ 


udies, the writer can make the following generalizations: the: 
a 


i 2d) For gate openings ranging from near zero (gate closed) up to 15 or 20 — 
o_o of full gate, with discharge coefficients related to the head on the 
_ weir crest, the application of the modified weir formula (first equation) re- _ 
ee in discharge coefficients S that vary from unrealistically high values a a 
zero opening to values nearly to the free discharge coefficient at 


For th the same range in openings, with — 
_ head on the gate sill, application of the orifice formula results ina nearly _ . 
ents. ‘This observation it is 


for given hee headwater condition, is a follows: tot 
Determine ‘the of water on the gate sill at tite 
- sidering this depth as “full gate.” For the majority of downstream settings imi 
ie on ogee | crests, it can be assumed that critical | depth will occur in the im- bf f 
mediate vicinity of the gate sill. Any corrections in depth or discharge, due | ¢aui 


to end contraction, velocity of approach, etc., should be made at this time. ae 
a) Compute and plot discharge versus gate opening from zero opening to a cI 


~ approximately 25 percent of “full gate” by means of of the the orifice equation and 


3) Draw a straight passing through the full gate point: of step (1) and °F 


tangent to the curve of step (2). This the curve for a 3 tigh 


conjunction with step (2) above, it should be stated that, although the mer 


definition of gate opening shown in Figure 6 represents good academic prac- [ alw, 
tice, a | good deal of manipulation is required for the precise determination of — era 

any gate opening. It is suggested that computations could be simplified 

oe _ expressing gate opening as a vertical distance above the gate sill, and by ad- | 
_ justing discharge coefficients accordingly; any differences in gate settings _ ade: 
and/or crest shape then would be accounted for by the approach angle 6, sus ; — 
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In reference to the author’s the results 
_ shown in Figure 7, it should be apparent that, for gate sills located down- 
stream from the crest, the approach angle consequently, the coefficient 
_ of discharge should be less than for a weir r crest setting. ‘This trend is amply _ 4 


| illustrated by the author’s Figure6. 


= ee The author’s dream of large models set up for | the sole purpose ¢ of re- 


os 


rms 


~ gearch in the field of gate discharge relationships: is heartily shared by the oe : 
: writer, whose interests in such matters » unfortunately, are not supported by , 
the means for carrying out such a program. . It should be of ir of interest to learn = 
_ whether the Bureau ‘contemplates | such a Comprehensive prog: program in basic 


‘The results of the pressure profile s studies are very well explained nf the a" 


Ag datum in the series of tests. In pansing, # it is of interest to note that the = oF 
spillway at American Falls, which was designed a good many years ago and i < 
_ without the benefits of present-day knowledge, is in no imminent danger of 

failure due to negative pressures on 

oa is obvious that this portion of the paper deals pr primarily with seals that at a 
are currently being developed o or r applied by the Bureau of Reclamation. It 


seal types and methods of application as there are gate designers. tre — 
merely describing several popular combinations would the 


| limitations placed on the length ofa paper, 


_ Although it is universally conceded that low head tainter alan should be es % 


> e equipped with side seals, it does not necessarily hold that they must also be pe 3 


_ provided with bottom seals attached to the gates themselves. In a good many a Gg 
early installations, the bottom edge of the skin plate was arranged to bear on 7 a J 7 
a creosoted timber, which was 8 bolted into the supporting structure s so that its &§ 


“vertical - lift spillway gates sealing against timber sills were installedinthe 
“early 1930’s; the original timbers continue to provide practically adrop- = &§ 
tight seal after more than 25 years of service. In a number of fairly recent Bay 4 
installations, the skin plate rests on a soft metal, such as babbitt or lead, ee ‘ 
which has been melted and poured into a structural channel, or other suitable — aS 


“always met with ‘the anticipated degree of success; their shortcomings gen-_ 


1) In cases concerning small gates, the weight of the gate » does not provide 


@ 

“adequate indentation of the metal, with the result that a uniform of 
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ser .. ihe Case Ol a rather deep gate, the hydrostatic pressure at the sill 4 ys 
| aclosed gate induces a wire-drawing action; this results in pinhole or shell a 
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rubber bottom seals are necessary, the writer anneniien die ‘not 


skin plate; he prefers to place | the seal on the e upstream face, so that oe aig a 
eof static pressure can assist in effecting final closure. In his opinion, the only — 

te practical difficulty that might arise from the latter practice is that an up- | 
- stream seal presents a local obstacle to streamline flow and to direct é 


In the application of side seals, the writer is strongly in agreement with iw 


the author in that a nominal amount of seal preset is necessary to assure 
a contact between the seal and the embedded seal plate. If no preset oq 


» 


1) Imaccuracies in field alignment of 8 seal 1 pl es and/< 


excessive out - -of - and sags. 


a my certain amount int of lateral pl play accompanies s the vertical movement of 


conse 
work: 


_tainter gates; it may be limited to a fixed amount by rollers or skids, ortoa | 


—= amount by limiting forces developed within the seals themselves Ly 


3) Field observations indicate that the contact face of a seal acquires some 


~ permanent set, resulting from age and usage; the softer the contact face ae? 
the more flexible the seal bulb (as for *J-seals”), the greater the amount ara 


Seal preset does, of course, frictional in addition to the 
; amourt attributable to hydrostatic pressure alone. This has led many to con- 
sider the possible application of retractible side seals, wherein all | 


an) forces: would be absent while the gate is being raised or lowered; full hydro- 


_ static pressure would be utilized to force the seals into contact with the seal ¢ 
a plates only after the gate has reached its cmanen position. Prominent es 


series of ‘model studies « devoted to ye ak the possible application of r re- es 4 
tractible seals to tainter crest gates. A significant observation derived from 
his studies is the fact that the retractible body of the seal must have con- Pies 

siderable inherent stiffness, so that the activating water pressure at lower 
levels can effect complete closure at headwater level, where hydrostatic uel 
a equals zero. - Moreover, the required stiffness would tend to Ree, . 
i es the flexibility along the contact face to an extent that could require extensive 

modification ir in order to regain a portion of the contact flexibility originally 


porates the simple “compression” mounting of the J-seal, in which the seal 
leaf is directly attached to the skin plate; it also takes advantage of the er 
greater leaf flexibility of the J-seal’s “deflection” mounting, while rect | 
venting the “outboard” mounting requirements of the latter. Perhaps the 
‘most outstanding advantages of the J-seal are its versatility and low cost; % 
several manufacturers of industrial rubber products offer a wide | 


ange in J-seal sizes and molded corners. The writer would be very much __ 


interested in learning how the L-seal negotiates the bottom corners of a = 
a His small contribution to an engineering handbook(9) has made the writer 


‘The Bureau’s angle seal (or “l.-seal”) appears to represent a compromise i 
between the two most popular types of “J-seal” mountings. . The L-seal incor-_ 
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Ba the limited audience would even pay the printer’s costs. rey 
‘The writer gratefully acknowledges the cooperation of Mr. Richard D. . 
Harza, A.M. ASCE, of the Harza Engineering Company, Chicago, , in 
consent to make public certain observations derived from the unpublished y es 
"works the F. Hi ; 


| 9. . Bowman, J. s., and 2. Pe: Gates” section of 
Handbook of Applied Hydraulics (C. V. Davis, editor - 
Edition; McGraw-Hill, New York, 1952. 
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“CALIBRATION OF 1 OF A SUBMERGED BROAD- CRESTED WEIR® ei. 


"Discussions: by P P . K. Kandaswamy, N. Rajaratnam, 


"presented a very interesting paper on broad-created weirs with upstream and 


downstream slopes. The new method of approach presented therein to find 
out the modularity limit based on the Froude — at the center of the — 


weir crestis commendable. | —— 
; all flow over a weir of any shape can be classified as modular or non- — 


a 


on the upstream head and is of the tail- “water head, its 


modular as explained below. When the discharge over the weir depends only 


the modular range, the over a broad- weir with 


cal upstream and downstream faces is given by the 


* B is the length of the weir in feet normal to the direction je flow and oe 7 


H is the head in: in feet o over the crest of the wei! weir. 


2, 3 ,4) by by various have shown that the coef- 
ficient C varies mainly with H/L where L is the width of the crest inthe aie 


of = is independent of the velocity of approach represented 4 


is the height of the weir crest above the Saal if 


‘The curve given in Figure 1: shows of C with 
H/L for a broad crested weir with vertical upstream and downstream Rae: poh 


iz When the broad crested weir is provided with upstream and downstream 


ae 


_ slopes, the flow pattern and hence the discharge characteristics are likely to 
_ be different. Considering the downstream slope first, it is known that when “ss 


water | flows over a broad crested weir, the critical depth occurs over the a 


crest. | ‘The e exact location of the critical depth section is difficult to define a as a 


can be seen from Figure 2. To be sure, the 1e stage of th the flow caheee ! from ? 
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ceases and the discharpe depends upon Upstream 
‘ i downstream head, its performance is said to be non-modular. The limit be- i ae a 
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Coefficient 


Woodburn L 2-00 Feet 
Besin L = 656 Feet o 


10 00 Feet 
hig 


g ‘the sub sub-critical al to 1 to the supercritical _ condition over the ‘crest. Hence, so long 
as the downstream slope provided is steeper than the critical slope for the err 

; _-Fange of discharges involved, as it happens to be in the case considered by ole 
the author, any change in the ‘super-cr -critical downstream slope is not likely 

to: affect the flow pattern over the crest. Consequently, the discharge coef- 


C in the equation CBH3/2 be independent of the down- 


"where is ‘known as the ‘ factor. ” The variation of y with H/L 
shown in 3 values of varying from 1/100 
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am 
il-water may bac be termed al 
involved, the tail which may be 
larity limit may scharge coefficient for swi 
the “modularity slope.” eneral way, 
1as been observed in a g int 
| has been o 
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Distances in Feet — 


__| Position Of Critical Depth 
Indicated OnEach Pra file 
Small 


TLL 

H. 


“ig: 2 WATER- ~ SURFACE PROFILES AND | POSITIONS a 


CRITICAL DEPTH ON WEIR WITH HORIZONTAL CREST 


‘to 1/1000. _ The results are taken from a series of experiments 3 conducted in in i ted 
thes Irrigation Research Station, Poondi on bye- -washes which are flush escapes > 


_ for floods provided in small tanks in the State of Madras. Further experi- 
4 _ ments have to be conducted to collect more information for steeper cael up 


> ‘The author has measured the depth of flow at a point located slightly — v 


responding Froude n number. ‘This has been used for defining the 
_ Limit, that is the stage at which flow changes from modular to non- -modular oy 
When conducting similar experiments at | this Research h Station, the 


-_ that sometimes the same depth occurs at more than two sections. “Hence, PT 
_ the writers feel that it may not be very reliable to base any s solid conclusion ae 
on the depth measured at the section mentioned previously, 4 
Considering the functioning of a broad crested weir in the modular range, i a 
any increase above the crest up to the modularity limit does not affect the rang a) 
upstream | level. This apparent anamoly has been explained by the writers(5) pag a 
toa the conception of the critical depth section occurring as a control section 7 
_ over the crest and of t the celerity of the breaking ‘surge. In the paper referred 


4 above, it has been concluded that until the modularity limit is. attained, the 
_ critical depth occurs somewhere on the crest. In fact, it can be seen that ae. 


_ witha a steady increase in the tailwater level, the critical depth | keeps on ow 
a “moving ng downstream over the crest up to the modularity limit when it falls 


a ‘directly in a qualitative manner with that of the author that the change in con= 


trol: occurs when the Froude 1 number is approximately 0.78 at the section 


<_< 
| 
| 
| 
| 
4 | 
4 
| 

1 
ot 


1220 er, 1959 


\ 


HE 


rested weir in the non-modular range — 
al 
reduction factor which is a function of the submergence ratio = and 
_ being the height of the tail water above the crest. we 

1e variation of the reduction factor f with the sub: 


sve 


Tres 


tee = 


is the modular discharge for the given upstream head and f isa = | 


4 
ah 
_ ir called Kodiveri Anicut(5) is given in Figure 4. Studies on aie ae on 


— 


= 


i 


it 


3 


where Su is the upstream slope studied is the downstream slope of 


The ante nature 


— 
— 
ram 


Engineers, Ist July 1952. 
King, H. i. W. ‘Hand-book of Hydraulics. 7 


«4, ‘Smith, D. «Open channel water measurement with the broad- crested 


Rao, J. V. Annual Report of the Irrigation Research 


; ‘literature on the broad-crested weir coefficients was incomplete. The 
analytical relationships worked out by the author between Froude number, 
degree of submergence and weir coefficients (illustrated graphically in 
* ae Fig. 3 and 4) give a sound basis: for selection on of | proper coefficients in parti- 

The writer has bess concerned with selection of weir coefficients for es- 
by _ timating discharge obtainable if reservoir dykes (rock fill or earth fill with | 

ig ns _ adequate top and slope protection) are designed as broad-crested weirs to 7 
v pass flood flows. From Fig. 4 it appears that for weirs with side slopes 1 fs 
on 1 or flatter, a coefficient "C’ = 3.0 may be used for submergence up to 

f he i 70% and where there is no likelihood of tail water rising above / weir crest, ‘ 
_ this value may be raised to 3.2. Perhaps it would be necessary to reduce ~ 

= re these values somewhat for design work to allow for variation of roughness" .s 
coefficients of the weir materials, velocity of approach etc. etc. Laboratory 


work done by L. W. Long and N. W. Conner(1) indicates that at high flow _ = . 


yates the discharge coefficient of broad-crested weirs is lower for paraffin i 
- surface than for glass. ‘The writer would be interested to know to what extent 


Smith’s values were modified for in the design office to 
7 4 ‘Figures 2 and 4 indicate that the effect of i increasing side ane of — 


crest paath 1000 feet and total head 10 feet, discharge capacity © = 134, 500 c. 


for 1 on 2 slopes and 127,000 c.f. s. for 1 on 1; % difference of discharge = = 5% ; 


Only). . This leads to the useful « derivation that flatter side slopes do not in- 


crease discharge rating significantly and considerable savings may be made 


& building weirs to ‘say 1- 1/2 ¢ on 1 1 slope i instead o! of 1 on 2, or r flatter, — 


Ss on Broad- Created Weirs with two different surface adhesion Condi- 


tions”; L. W. Long, N. W. Conner. North 


to 


Engr. +» The ‘Shawinigan Eng. Montreal, Canada. 
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TA TABULAR SOLUTION OF OPEN FLOW LOW EQUATIONS 


, R. Ivan M. Nelidov, 


SILVESTER .—The has provided a very useful table giving 
7 values of specific energy readily and accurately. The author’s equation 3 has a 
appeared previously. (1) The analysis of the dimensionless form has been 
restricted to the rectangular can be modified to o include other 


where by is the water surface ~ and “st w is the | width h of a rectangular 


channel equivalent in n area to th to the c values 


for trapezoidal channels for any operating. ‘For 
shapes a scale drawing could provide values close enough for most practical — 


problems. Where depths are being determined an initial estimate based on ~ 


the original equation could be the basis of a correction applied for the shape 


j aa This procedure may may appear a little arduous in order to obtain one or two 
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Csr 


= 


where dx is the lenath ht is friction 


. ‘e specific energy can be determined from the table supplied by the idee, 
rectangular channel this is done ‘directly. For other shapes the 


li re The correction in fact is only applied to the velocity head term and, as can f 

be seen in the author’s Figure 1, this is small for mild slopes once d/d, > 1.5. | 
te - Below this value down through the critical depth and for steep slopes the -~ 
rection assumes importance. The velocity head is found by subtracting 


= the value E as read in the table. Ar An | example should assist to 


in channel of bottom width 25 ft, side s slopes 10 on 1 _— and bed slope "aed 
A 0. 00079 delivers 900 cusecs into a reservoir whose depth at the junctionis 


. ~ = 7 It is required to find the distance upstream at which the depth is within 


percent of normal; Manning’s n = 0. 025. 


_ Graphs nomograph for ease il in normal and criti 


de 


4 . Thus the total length of the backwater profile is 3445 x 3.28 = 11,330 ‘ft. “i 

_ It may be remarked here that the same problem was solved by the conven- 3 
a _ tional step method in reference 3 where the answer given was 14,360 ft; but a 

_ it were reworked using the average values of friction and kinetic energy i 4 
terms for each channel ‘section, as above, the answer would become 11,300 ft: 

~ Actually using the table without the ‘shape factor correction results in a length 


= 


o 


a. 
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IVAN NV M. NELIDOV, ! F. ASCE.—The writer was attracted by the method 
esented by the author for facilitating solution of open channel problems hi aol 
In accordance with the preference of the writer for graphical solution of 
draulic problems he wishes to a of this 


oblem for trapezoidal channels. 


Since for critical flow the a fa head is h 
2T. 


pr 


curves es representing of and B/m are plotted ¢ on Fig. 2. 


‘For ‘the specific fic energy 1 to critical the expression 
can b be 


2s| 1.0 | 2098 2.000 | o0¢ — 
— 
~| 
9.0 | 274 | 2807 | 0.067 [050 /0.038 |2.778 ot 
~*~, 
— 
| 
) ft: 


VE 
AS, | 
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where: 


+ 


representing 


and ‘B/md. The application of these c curves will be illustrated on the. examples 


ae Example 1. Data as er author, except that the side slopes a are 1.5 to 1, or a 
m= 1.5. Then, Qo/ g1/2 = 400./5.67 = 70.7 and B/m = 9./1.5 = 6. 
curves ‘on Fig. 2 we obtain d, = 3.25 ft. and E 25 = 2. 15. With 

obtain d and = 2.11 curves ‘on Fig. 3. This 


us d; = 57 x 3. 25 = 1.85 ft. and = 2.11 x 3. 25 = = 6. 85 ft. 
Example 2 2. Data as per er author, except side ye factor m = 2. ‘Then 
| 1/2 = 300. /5.87 67 = 53. and 1 B/m = 10. = = 5. From curves on Fig. 2we 
obtain = 2. ft. and = - 2.4. With B/md cl= 10. 0./2x 2. 5=2. 


ao) 2.43 and Eg) = 2.4325 = 6.07 ft. 
For the 1 nar 


row section, wi! with the side slopes, we obtain: with 
88. from. Fig. 2 and B/m = 9/2 = = 4. 5, dog = =2.65 ft.; neglecting losses, “ 
5 


Fata 6. 07/2. 65 = 2. 2.29. Then, from Fig. 3 and with B/md, c2 = cecal 6 


1.7, we obtain = 2.32 and dy = 2. 32x 2.65 = 6.15 ft. 
Example 3. data by the author and factor m= 5, also with 


= 5. 5./2. 95 = 1. From Fig. 3 with = 8. ./BX2. 95 - = 5.43 we find 


=1. 1.83. =1. 5/2. 95 =1. 83 -. T= 1.66; with th this we 
1.48 and dy = 1.48x 2.95 = 4.5 39 ‘on 
circular conduits if we consider the central. angle a a subtending the 
water width T, the same will have an expression T = DSin 1/2 a, where D i. 
* _ the diameter of the conduit; the » depth of water will be d = D/2(1-Cos 1/2 a) 
and the area of cross section A= D2/8 (a - Sin a). The U. S. Bureau a yy { 
Reclamation tables contain coefficients necessary for plotting curves similar 
_ to those for trapezoidal channels. (Table 12a, p- 48a, Ed. 1923). The rela- 
a tionship between the critical depth ratio to diameter d./D is given by - 
4 ficients Q,/D°/2, which we will ll £3(d,/D), so that Qc = fi (d,/D)D*/2; 


4 relationship between and the area. is by the coefficients 
48), which we will call f2(4/D), 80 
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‘multiplying by d, and substituting 
‘This cubic always | has two positive roots | or depths (that is, 


when Es “and the constant term are positive. A nomographic solu- > vi 


tion for the enteoction of these two positive roots is mien on — we? 


ine UZ 


WY 


—! ation which arises in 
hydraulic probiems Can be Obtained Irom the author's equation 
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the diagram c the -author’s 
0 ple. The equation to be one is 
nce the scale of 4 


A straight edge through 3 
curves Es 3 
= 


ls 


that the pertinent cubic equation is. 4 


+ 


_d 


This chart ¢ can be constructed to various 
and desired and the 5— scale can calil te 


5 a the nomographic solution of the ne hydraulic tien 


: 
‘The 

 eithe 
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— 
— 
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dimensionless specific-energy function, in either graphical or tabular form, an L 
can be a useful tool in the solution of certain open channel flow problems. — oe. 
The fact that the method presented permits the inclusion of energy losses, , a 
f either as a numerical value or in terms of the velocity head, gives ita ae 
& tinct advantage over previously published procedures. This was pointed out - 
> to the writer by the author some time ago when the latter was a graduate 
student in a course in open channel flow given by the since that 
_ time the method has been included in the course material. — ie “sa 
_ ‘The use of the diagram and table of the paper are a to the rec- 
wae cross section, and the method cannot be used for the trapezoidal a 
= because the ratios involved change with the geometry of the channel. — 
There are, , however, two other shapes, less commonly used, for which a simi- 
ilar procedure may be used. These are the triangular and parabolic channel 
| cross sections. Since the determination of the depth for a given energy for 
¥ : b these sections is even more difficult than for the rectangular channel, con- 
t siderable work can be saved through the use of the author’s methods if — 
ae general criterion for critical flow is found in most textbooks and is de- 
rived by the substitution of Q/A-V in the author’s Equation 1 and setting the 
: derivative of the expression equal to zero for minimum energy. _ Noting th that 
— the width at the fr free surface, results in the expression ane 


"which is satisfied in any open channel when the depth iscriticak | 
@ Consider now a symmetrical triangular channel, vertex down, ¥ with a ver- 


: tex angle 8. Both A and Tin Equation 9 men be veaiiaiein int terms of of d and 


sr 


? 


The equation for specific-energy is 


af 
can by the same procedure to be valid an u 


The channel may be defined by the equation Tekal/2 2 


bait 
| 
“4 
| 
— 
| This is reduced to dimensionless form by following the author’s procedure — 


and the dimensionless apecific- 2 


is evident ae the section, either ‘Equation 12. or 
; _ Equation 14 may be used in graphical or tabular form for the solution ae 
4 ‘problems involving their respective shapes and that the curves ‘which they | 
represent are independent of the values of 8 and 


ae may be of interest to note that the pressure- momentum function, : fre- = 
a, used in connection with the hydraulic jump, aay also be expressed in 


and making the substitution Q2 = a3 b2 g ag and dividing by d. gives the 


: 


channel is 
_P+M 
and for the sunita section 
= 


Ta 2 
— 2 
4 
> 
— tal 
| 
— 
: 
— = 
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— other channel being similar in form. 
eS - Knowing either of the depths, say do, for ahydraulic | 
Ge ratte is calculated. The ratio is then selec = 
having the same abscissa value as do/d, and the value 
— Alth ugh the conjugate depths for a rec angular 
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xe 
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rt 


igure Dimensionless pressure-momentun 
may be other means, this is not true for the 
‘ shapes. The conjugate depths for the latter can, however, be easily deter - t 
- mined from a graph such as Figure A drawn to an appropriate scale, or: from 
table the at a higher degree of accuracy. HH 
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Discussions 3 by K. Henry J. Miles, 
MB. McPherson, Nicholas Bilonok, and J. V. Dawsey __ 
n, 


= szESZTAY!.— —Recent investigations into the movement of fluids p pr 


‘ “com common theoretical basis for the analysis of flow in closed ci conduits and in 


{ duits are compared in a very clear manner to the inclusive formula and all 


gram of Colebrook-White, , respectively Rouse-Moody in the paper of 

_ The roughness of open beds is usually expressed by the n value of Kutter 
- Ganguillet, while the equivalent particle diameters k of Nikuradse are used 
in theoretical investigations and in design computations of closed conduits. | 
Combining the Karman-Prandtl formula for rough pipes and the Manning 
formula H. Rouse developed for the relationship | between the two ‘roughness | ey 


ool . : : « In order | to study the effects c of fact factors which cannot be taken into account 


analytically (kinematic viscosity vy and the gradient S) and errors in observa- 


: ¥ tion, a statistical investigation has been carried out recently. ‘The evaluation 
. e of the 116 data obtained experimentally, respectively by measurement, and az 


extending over the entire turbulent range in both open beds and closed = 2 
duits, yielded the of two variables (Fig. 1): (3) 
7" 


mentioned above, the statistical investigation f failed to reveal the effective of 


cist 


a In his paper ‘Mr. W. Moore points toa very ‘important prerequisite for the 


widespread a application of the developed fi formulae, which 


a. Proc. Paper 1962, March, 1959, by W. we Moore. 


he, 
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4 4 mations can be eliminated. An easy solution can generally be found to ee 
_ problem in case of conduits having a circular cross section, all characteris- 
tics thereof being the function of a single parameter, the diameter. In the 
of the most common type of open canals, ) the cross 


Graphs were recently issued by the Institute for the ulic 
_ design of canals and pipelines. (4) _ Attempts have been made therein to give a 
_a direct solution for the dimensioning of the optimum trapezoidal canal cross” 
section for any given Q discharge, terrain and 
The problem will not be unambigous, unless certain characteristics 
: “ the trapezoidal cross section are either assumed as constant, or 4 
sidered as a direct function of the water depth (or any other dimension). ae 
= basic assumptions for the construction are shown in Fig. 2. ies set Boe q 
_ Variations i in the shape characteristic _B = a/h with the waterdepth h, have 
i wi adopted after studies on stable beds (Lacey, Lindley, Girskan and others) 
_ and on basis of empirical data. The relationship between the bank slope 9 g 
and the waterdepth h represents a statistical summary | of solutions common- 
ly used for unlined beds. _It should be noted, that the method to be presented 
_ below can be extended by the construction of an adequate number of graphs to 
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le is a the ¥ wetted perimeter, A ts on sectional area), aunndine the 
a variation of the hydraulic parameters of the cross section with waterdepth. 
th the a expressing the conveying capacity of me beds 


32 10g (00 0676 
ers) 981 the constant of gravity) asa transformation 
Colebrook-White formula by E. S. Crump. 
The inclusion of the second term in Eq. (5), expressing the effect of vis- ia 4 
sa cosity and slope, is not necessary in case of conduits and canals carrying te i 
water unless k < 1 cm, respectively n <0,017.(3) 
3 Consequently, with the substitutions ; according to Eqs. (3-d), (3 -e) and (5) 
the by E (4) is - for given roughness — 


= 
~ 
4 ‘Having sel of B and 9, respectively the corresponding fil 
condition equations, the Construction of the graphs constitutes a 
_ | hydraulic problem. The first step is to establish the relationships eo 
— 
) 
7 re 


tively k = 7 cm) is shown in Fig. 3. When using the graph, the mean flow Ae 

velocity v and the waterdepth h is determined for the given discharge Q 
and surface slope s, whereafter the B ei values corresponding to the ed 

waterdepth h can be obtained from Fig. 2. With these values available, all 4 


‘es 7 further parameters « of the cross section can be computed by Eqs. (3), or can ; 
> 


obtained from the corresponding graphs. 
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: aluce obtained from the graphs may be used as a first approximation — 


o en i if local conditions asr regards roughness, bank slope, or the the ratio ya/h 


__ By compiling a series consisting of an adequate number of sheets, the - 


fe solution according to Figs. 2 and 3 can be extended for any c combinatior 


Wiley Wiley and sons, New York, 1980, p. 115. : 


2 Crump E. S.: : Flow of fluids in conduits and open channels, Proc 
Aut 1956 pp. 
Szeszta K.: Az aramlasi sebesseg szamitasa (Computation of on. J 
velocity) Viziigyi Kézlemeények (Hydraulic Engineering) 1957. No. 1-2. 
(in Hungarian with summary in English, French and Russian). 


4, Szesztay K. Tervezési segédletek csatornak | és csévezetékek 


_ méretézéséhez (Graphical aids for computation of flow in open channels = 
and conduits) Vizgazdalkoddsi Tudomdnyos Kutaté Intezet, Tanulmanyok | és 
Kutdtasi Eredmeények (Proc. of the Research Institute for Resources) 


HENRY J. J. MILES, LES, 1 | M. ASCE. .—Professor Moore’s paper isa 


vil contribution to a better understanding of the subject of pipe friction. The — a 


_ writer is particularly interested in this paper as he had independently de - 


rived a similar relationship between exponential formulas and the general Hs, 


_ resistance diagram, and for many years in his classes in hydraulics has been 4 


ia using this idea for introducing the overall concept of fluid friction. As the ih 
T ~- points out, the large number of exponential formulas is not only con- © 


fusing to the ‘beginner but also gives the wrong impression of the inter- a: 


“cae sn of hg, D, Q and V unless it is realized that each formula nomi 


only to a limited range of Reynolds number. 
writer used a slight | modification of th s formula 


‘In in the range of laminar { plot of 


‘general resistance diagram with a negative slope of 1. Therefore, m = 


or high Reynolds numbers the slopes of th of the lines approach O, m approaches — 


Prof. of Ci iv. » The University of Arizona, ‘Tucson, Ariz. 
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September, 1959 


the viscosity term having disappeared completely. inter mediate 
umbers, as the author points out, the value of m varies from about 1.7 to 4 


Reynolds number increases: hg varying initially as V and inversely as i 
nae and gradually « changing to the relationship of varying directly as V2 and in 
_ versely as D; viscosity initially playing an important role but gradually de- FF 
= _ creasing in importance until it finally disappears from the equation. It is also 
apparent from observing the slopes. of the lines on the general resistance 
y _ diagram that the rougher pipes approach complete turbulence at much lower *# 
Reynolds numbers than the smoother ones. 
‘ ae, It will be noted in the writer’s formula (1) that the sum of the exponents of — 
D and V always totals 3. This tendency can be observed in the various ex- 
es ponential formulas yumpesied in the past by Hazen-Williams, Scobey, Barnes, 4 
Lea, and others. Since the sum of the exponents of D and Vv in these ex- 
® aay ponential formulas does not always total exactly 3, the question arises as to aa 
or whether the discrepancy is caused by slight errors in the empirical formulas, 
whether the idea of { representing corresponding values of f and Np on 


In conclusion, it is interesting to use the methods proposed by the author to’ 
compute the numerical r relationship the woe ratios 


yer of D are adjusted slightly so that they total 3, the 
quated to the writer’ 's formula (1) 
Sar 
from which it is s obvious sia mn must equal 1.84. If this value of m and the 


kinetic viscosity of water at F. are inserted in it reduces to” 


& 


sine 


— 
0.16 
phe 
| 
a 
4 
— 
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it corresponding ' — of f and Nr are on from the § general resistance 


for each e/D line at the point where its. negative ‘slope is s 2-m 

0. 16, and substituted in this last equation, Table IV is obtained. if these po 
"values are plotted on semi- -logarithmic coordinates, the resulting curve eap- 


for Water at 60°F. 
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_M. B. McPHERSON, 

sentation is directed towards a a interpolation « of the resistance 
diagram,’ it is necessary to review timmetions and assumptions which must By 
be considered in using the ram. : 


~ Colebrook, ( 5) ) noting that the general trend for new commercial pipe oe 


was a continuous transition between the characteristics of hydraulically os 

I Water » Philadelphi 


— 
— 
is 
— 
— 
ag 
ia 


“smooth” and artificially “rough” pipes, arrived at an asyupietiec “transition 


q ae law” of good approximation based upon integration limits for these extremes 
: a given by Nikuradse. The Colebrook “transition law” is graphed in Fig. D-1, 
a photo-copy of Fig. 1, Pp. 136 of reference (5). 
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in which 5% (much larger size than the remainder) was cemented in a diamond 
eer pattern | in a test pipe by Colebrook. _ The original raw data for the three a 
ee os classes of new commercial pipe were selected from a variety of published BS 
_ American and private British references. Colebrook determined the equiva- — 7 
~ lent sand grain roughness, e, for each individual commercial test pipe direct- 
ly. and indirectly by means of the Karman-Prandtl “smooth” and “rough” pipe 
relations. From these calculations for each pipe, he then determined the oe 
“mean” e for each class. _ The extremes and ‘Mean for each class are sum- 


roughness. has a variety of test data (mostly 
= diameters) with the Colebrook “transition law” » in his Fig. 110. It is to be 


noted that Colebrook interpreted his results in terms of a system of transi- 
tion. curves, whereas only his semi-empirical “transition law” has been in- 
corporated in the general resistance diagrams presented by Rouse(6) and 
4 ie ae as being “a close approximation to the actual resistance law. »(6) bi 
_ Inasmuch as the majority of the tests supporting the validity of the transition J 
- region curves of the general resistance diagram were for new pipes of small — 
diameter, there remains some as to the probable accuracy in the 
we of this diagram for large pipes. 


‘Ries 
‘For the Types Considered by ook: 
Galv. Iron” 0.006~in.(.0005-ft.) 
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‘The mean G. I. transition curve of Fig. 1 would lead to values of f up 
oe to 5% higher than the “transition law.” If the range of calculated e values” . 
i be) for G. I. pipe given by Colebrook (Table D-1) are applied to the general BR . 
e sistance diagram, Figure 1, the variation in f about that for the mean e is _ 
i, of the order of plus and minus 15% for a 4-inch pipe at a Reynolds number of :) 
105, This isa reminder that the first four equivalent uniform sand grain 
a, roughnesses for new commercial pipe tabulated in Fig. 1 are meanvalues. _ rs 
: * (Ranges of values, such as those given by Streeter, (18) would seem to be less 
_ “4 misleading). "Most references suggest that the accuracy of estimating f ae 
the “rough” region (for which there is a paucity of supporting data) “rr . 
_ probably within | about plus or minus 10%: for “smooth” pipes the probable > be 
variation is about plus or minus 5%. Recall also that the supporting data for 
ae _ the transition region analysis was selected from laboratory and field data 
considered representative of new, clean straight pipe without fittings 
‘ (Colebrook made arbitrary head loss deductions for bends and fittings from — 
_ The universality of the “smooth” pipe characteristic for ‘steady near- ieee 
incompressible conditions has been demonstrated using oil, air and water 
_ (reference (3), p. 203). Nikuradse’s experiments as well as all the com- > 
_ mercial pipe tests used by Colebrook and subsequent adaptations were con- | 
: - fined to water flow. The trends of the functions for new commercial pipe with 
» 8 random or ‘statistical’ type of roughness are indeed remarkably similar. ak ; 
“ be ‘However, the preciseness of the general resistance diagram implied by the 4 
: an “exponential equations” of the author does not appear to be justified. The 
ie author leaves the impression that the characteristics of new commercial nine 
are fully established and completely understood. The extensive studies = 4 


artificial random sand grain roughness by Colebrook and White(28) and 
Square-groove machined-roughness patterns by Streeter(19) indicate that 
more remains to be evaluated. -Rouse(20 


warned some time ago “Sy 


describe a any roughness other than uniform. In reference Rouse recom- 


mended the Colebrook relation “as a close approximation to the actual re- : 
sistance law.” 


_ states that the effect of commercial pipe roughness “has not yet been ex- 
Blained in terms of our of the action of artificial ” The 


writer | assumes that the in n fluid tendered 
“American undergraduate civil engineers for more than a | decade has left at aa 


- least a residual impression, and is confident that contemporary designers i 


limitations inherent in commonly used empirical relations. He assumes 
Te that a comparison between the Hazen-Williams formula and the general re- ee ; 

_ sistance diagram, similar to the presentation by Vennard, (22) is offeredto 


21) 
.” From the standpoint of velocity distribution, Robertson(2: ) Pi 7. 
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remarks are not intended to disparage the vast and invaluable knowledge ac- 
quired to date, but to warn against the false sense of exactness implied by the h 
too-literal interpretation of the trend-curves as in the author’s “exponen- 
answ 
&§ _ twen 
abou 


"mechanics in this ‘country. There has certainly been a visible 
_ grading in American design practice since that time and it is hoped that the e.. ‘ 
EE author’s pessimism does not truly reflect ‘prevailing conditions ¢ except in the | * 
extreme. The relatively contemporaneous heated discussions on formulae to — 
_ which the author alludes(10,14) are specific to British practice; some of the fe 
in these references are peculiar to British usage. 


usually a appraise » the | expected initial, or “new,” 
a of this type of main by means of the general resistance diagram. i 
Since designs are normally based upon ‘future’ requirements, an estimate of 
both new and aged roughnesses is necessary. - Recourse is made to well- — - 

- documented and carefully selected data such as presented by the U.S.B.R.(23) _ 

and/or t to local re records in or critical (design-wise) 
author has made some valuable and. timely comments the 
ing “Equivalent Pipe Concept.” However, the numerical values of Table 

used to support some of his arguments, pass beyond the realm of design prac- = 
_ ticality. In water works design, an overall arterial, feeder, or express main 
head loss exceeding about twenty feet per mile (equivalent to nearly an 

j _ pound per city block) at peak flow can seldom be tolerated. _ ‘This criteria is 
| equivalent to a flow of less than 3 cfs in the 12-inch pipe and less than 0.5 cfs 
“a tn the 6-inch pipe of Table II. Even ignoring the unrealistically excessive ry 
4 te losses cited, it is implied that the field data from which the pipes would be 7 
E would be sufficiently precise to assign three significant figure = | ‘ 


4 


coefficients and exponents to the “exponential equation.” Such precision is a 
_ rarely attained and should not be expected for typical distribution a 
a's Equation (20) represents nothing more than a restatement of the equiva- “<= 
lent pipe concept and has been used for a variety of purposes in this or simi- a d 
lar form for several decades. The equivalent pipe procedure is occasionally f- 4 
a employed in modern network sualeen as a computation aid (for instance, to es 
reduce two pipes in one street to a single resistance for simplification, if <4: i 
A justified). , Subsequently it will be shown that the procedures proposed under — 
the heading “Exponential Equations and Complex Pipe Systems” have little PI 
P place in the analysis of distribution networks. _ Further, the writer submits a 


en _ that when a distribution problem is as simple as the “complex” cases fs 
cs _ vaguely mentioned by the author, there is no logical excuse for using anem-_.. 

a aon fit to the general resistance diagram except possibly for a first trial; q 
a to be consistent, the diagram itself should be used in determining the final E 4 


- answer. The author states that with a twenty-fold variation in Q (quite hy- ot a ‘ 

% pothetical) the exponential equation will agree with the diagram within a am 

| For his first example in the Appendix, this statement might appear to be cor- 4 q 

| rect: the head loss for the 0. 544 coefficient given would be about 0-2% below 7 
; A that for the diagram. However, there is an arithmetical error in the example — ae 

| and the head loss should be 4. 88 -feet and the coefficient should be 0.454. A Be 3 

| twenty-fo1a increase in Q for the corrected equation will yieldaheadloss =~ 

about 17% below that for the diagram. The Cc for the corrected example a 

145, which is not likely to be associated with a small, 12-inch main. 
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Regardless, a more comparison would for a doubled flow . 


_ which the the corrected equation would lead to a loss about 2% below that for the 7 
diagram y while aC of 145 in the Hazen- -Williams formula indicates a loss 


s ‘periods of near- -average demand a and since peak hour flows in large srterial 
_ systems are typically two to three times the annual average, a rate approxi- — 
mately double the test rate should represent a estimate of the 


a Distribution N Network Design 


Network over 90% of waterworks distribution 


a — design problems. At the present time there are two rather widely used me- a 


chanical means for the detailed calculation of head losses in complex gas and and — 


water distribution networks: by direct analogy using the McIlroy Network ob 


a ie. Analyzer, ( (24) and b iteration \ via the ‘Hardy Cross relaxation procedure using 
a al i. Calculations on a digital c computer are restricted to a constant ‘m-value 
(1 85 normally used for water service), since a large number of exponential — 


equations 2 are being solved ‘simultaneously. Introduction of of a varied relaxa- 4 


tion factor (differing values of m) would hopelessly prolong the iteration ease 

_ process. With the McIlroy Analyzer, the m for individual pipes can be ar- sat 
: to fall | anywhere between 1. L.85 and 2. - Thus, in an analogue | calculation, 
it would not be impossible to as assign different values of m to individual pipes. 

ee a _ However, the purpose of these calculations is normally for the designof _ 
- facilities to meet peak de; demands under future conditions (head losses for peak 

oy demands under existing ; conditions can be estimated more easily , and more “en” 
"expeditiously by means of extrapolation of direct field tests). In this connec - 


ion it would p perhaps prove fruitful if, some e of the nor 


J 
ig "Neither allt the arterial pipes ina system nor rT segments of the pipes 


‘calibrated’ are accessable for field head loss measurements. Field head 
| 
loss measurements can seldom be made over a sufficient range of flow, even 
-- induced, to cover peak demands. (Design of large arterial networksis 
generally governed by peak-hour conditions). Field head loss 
: over a range of significantly different rates of flow are very expensive; _— 
_ measured values normally would have to be extrapolated to greater rates for c 


bp, After ‘skeletonizing’ the system to an arterial network, loads (or ‘draw- 


iis location. The . proportion of the total demand as assigned each load is based od at 


upon annual averages. (The most reliable information, meter readings and - 
leak-and-waste- -survey data, i is usually in in terms an demand and rate). 


Due to the absence of adequate information, itis commonly 


| "sumed that each local load fluctuation has been and will be proportional to the 
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= = __ cation, proportion and variations of loads for the future are all necessarily = ——s 


and additional but less important approximations would thoroughly 

discourage a a designer were it not that ‘control’ analyses of existing networks 
_ performed under these assumptions usually check out remarkably well with 

_ field system-head-loss measurements. One aspect of this fortunate coinci- 

7 dence has been at t least partially explained | by Mcllroy. (26) He made a a 
of the effect of random errors in assumptions of values of C orf for the odes : 

_ individual pipelines in networks and found a “pronounced tendency for rather hes 
Hepes randomly distributed, individual errors to cancel themselves in the 4 
_ determination of total head losses.” ’ He conc concluded from his study that ‘ “there 

appears to be no strong reason to insist on the use of rigorous expressions in pen | 

= solution | which is not based on precise data. ” And nd further, “For network ~ lial 
analysis work, the widely used Hazen-Williams formula is entirely y adequate, .- 


if the same Samee of wane is ee to the ae as is necessary in the 


; tele the writer is forced to reject as obviously unrealistic (and clearly not 
_ proven) the statement by the author: | “Although the methods for analysis of | 
networks will not be considered here, it is pertinent to note that febceet 
resistance relations will be most convenient for use in this sin of f analysis.” 

The Hazen-Williams Formul 


sideration in the of express or feeder a 
iE m-value is not only convenient but practically necessary for network 


— 


analysis. ‘This equation can be compared t to any point on the general resistance a. 9 
in the turbulent flow r region. The statement “Thus the Hazen- “Williams 
formula with its exponent of 1.85 corresponds to a combination of relative ca ari 9 4 
roughness and Reynolds number lying along the line m = 1.85 on the Moody ae: a ; 
diagram” is completely incorrect. - ‘Fig. D-2, a photo- -copy | of Fig. 5, p. 258 
of the author’s reference (14) clearly shows ‘the comparative relation between ae 7 q 
and The equivalent friction coefficient for the Hazen- formula 


ior yates at 60° F. Comparing Fig. 1 and Fig. can seen en that the 
i. curve in Fig. lis fora variable ranging t about Cc = 140 


_ The author’s Fig. 2 is a redraft of Lamont’s Fig. - y The table which enld 
have accompanied the author’s Fig. 2 is given here as Table D-2. Lamont did 
‘not give the original data. All but twenty of the ge toa cited have been ‘a 
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on do not collectively describe a single-valued exponent. Obviously, neither can 
the Hazen-Williams formula be described by a single line on the diagram. 

Re in Fig. 2 happen to be for new G. I. pipes with a diameter ranging between 4 Sapa 
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Reproduction of "Table Summary of Records 1 for New Pipes", Pe 


of | Diam: oedty: Recom= 


: 0.350 02556 3,40x102 0.0016 
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0.24 0.000 
10.12 4.093305 0.0050 
46,00 28x10° 
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6. 35 and 2-inch. These curves, as 5 well as several other : irrelevant curves = | 
_ for diameters up ‘to about six inches, should have been omitted in the redraft = i: = 


__ by the author; effects arising from differences in manufacture on roughness | a 7 por 

 arel not considered to be separable for larger sizes, but small sizes, parti- 

cularly under 3-inches, are known to exhibit dissimilar ‘characteristics. In the 

_ terms of distribution design, any size smaller than 6-inches would be of little a 

interest. . (The curve well below w the “s smooth” c curve is for an 11-inch “ spun 4 _¥ : 

pipe). the curves for the small sizes omitted, itis pos- 

‘sibly reasonably to argue that “for (selected) new pipes as exponent m from | &@ 

80 to 1.90 is common ---.” In Colebrook’s(9) Fig. 2 (reproduced by the 

 Hinds(12) as his Fig. 8) which ‘summarizes selected tests on sixteen new 8 ar 

“centrifugally spun concrete” lined C. I. pipes 4-inch to 60-inch in diameter To 

_ and on six new “centrifugally spun bitumen” lined pipes 11-inch to 46-inch ies | _ = 

- diameter, the divergence from a “smooth” pipe characteristic is slight. The tel 

same result is indicated in Table D-2 for spun linings (but it appears that = bi ble 
much of the same data was probably reused by Lamont) jg .. 

i; a: - Although the selected data discussed above is ‘reassuring, it is not neces- i 

~~ sarily representative of typical field conditions. In Table D-3 is presented a 4 

Es summary of results from field measurements on cleaned existing pipe with — 

‘new cement linings obtained under the direction of Mr. Herbert Thamer, | wy 

os resident engineer in Philadelphia for the Pitometer Associates. Not manyof | P 

of 
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=: of a machine- applied lining. The scatter can be attributed to differences 4 
is local losses due to alignment of pipes and bends and fittings, assuming that 
tests were of consistent accuracy. 

a The author’s Fig. 3 is a redraft of Lamont’s ‘Fig. 4. (14) - Colebrook and 
-White(4) in 1937 had observed that “old pipes - - - - obey the square law a 
and that their graphs included rather flat curves in the transition region 


necessarily results to be from cleaning and a 
To plot “old” pipe friction factors on a general resistance diagram is to sug- oa A 
gest that the increase in roughness with age indigenous to a} given water | —  - 
tem (i.e. incrustation due to a variety of local causes) is as generally = 
3 ble as the somewhat inexact case of new commercial pipe. It is difficult to 
_ believe that incrustations are necessarily of the same pattern and form, or 
- yield equivalent friction results from one system to another for chemically © 
4 different waters. - Details on the data plotted in Fig. 3 are lacking in Bie “a 
; reference (14). The evidence presented by the author does not provide the 
— reader with an adequate basis for concluding with any assurance that “for old 
| pipes an exponent of nearly 2.0 is common. ‘4 A partial basis for this argu- 
‘ 3 ment is given by the author by means of the m-values of his Table Ill, in 
_ which the Hazen-Williams C turns out'to be only 89 for a 3-year old, and 100 
a 5-year old cement-lined 24-inch pipe. In both instances it would | appear 
that the lining work was ‘substandard, or a nuniber of bends and fittings were — 
= __ included, or that tuberculation in Austin is extremely rapid, or that the meas- 4 
" —— were not very accurate, or some ne combination of these factors. For 
data from the Pitometer Company for some “old” cement -lined 
e pipes in Philadelphia are presented in Table D-4. It is hoped that the a “<n = 


supply similar details the pipes listed in Table III. 


Consideration of detailed refinements for only one factor, pipe 
may create the illusion that aes default the other variables are associated ro 
 inings in a network, the Hazen- Williams formula is of the same order of re- - 
as as the general resistance diagram when all other factors affecting 
the design a are considered. ‘The convenience in network analysis s afforded by re 
2 a relation with a single exponent (relaxation procedures are made feasible; — re 
tables and other computational aids are available) outweighs the probable 
greater accuracy of the diagram for | existing | pipes in good condition which 
Fa have been field calibrated at one or more rates of flow. When “old” pipes — a q 


7 _ comprise | part of an existing network, the expediency of a constant m- m-value ¥ 
4 for flows in all p pipes appears to be justified since arterial mains in }poor caret 
must often be cleaned and lined to meet expanding demands for the 


| "an important at amount of flow and hence it’s contribution to the 8 system head r 


these tests represent a near-" Tsuc. There is considerable 
org ar in the vwolued of deagnite noscih consistent ‘statistical’ 
& 
ie a Continued use of the Hazen-Williams formula in the analysis of water dis- _ % a 
| __tribution networks will persist until such time as the accuracy in appraising a Se 
a 
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loss will be tenn minor. At the present _ tale is insufficient evi- 


- dence available to warrant adoption | of further design complications eo 
use of a mixed set of m-values for analysis when a network includes =, 
In the of non-network vk express,’ ‘pumping,’ ‘feeder,’ etc.) the nie 


DOr 
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_ one value for a given pipe is ‘mode to coincide with the diagram. . The author _ 
treats only the case where a precise value of e is known in a rather 
___ It is quite evident that there is much remaining to be investigated from the 
standpoint of theory and basic ‘experiment if a higher degree of precision in > ia 
a calculations is to be achieved. Periodic shuffling of warmed-over data i 


Conduits will s stimulate interest in farther research on Pipe, 
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tself. The author dis i 
sary substitute for the general resistance di a system of infinite ap- 
| 3 dains the use of empirical formulae, yet propos y oat 
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- (See Proc. Paper 1962 by Moore). ee 
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; authors general pipe resistance diagram | shows the effect of the various fac 


ong 
practicing engineers due to its inconvenience as compared with similar oo. 
exponential formulas for water flow where viscosity has little « effect and “ae 
_ restricted range of Reynolds number, | and the velocity of flow, also the dia- 
_ gram in its common form, trial-and-error solution is required for finding = 


; the rate of flow or pipe diameter; and for given boundary material the pipe 


| diameter must be known before | the relative roughness may be evaluated. nee 

| The diagram is based on the Prandtl-Karman formula for pipe resistance 

; -Telations which are improveu by the experiments of Nikuradse, Colebrook, 


iii 
> 
a 
— 
= 
—— 
se | neers. Valuable experimental data regarding resistance in existing water — — 
=— 
— 


r, 1959 


‘and White in form pines Table aise 


f Tred in Equation (1) is a combination of variables, 
_ which are evidently applicable to the problem, in very complicated etore 
7) ny _ where f is dependable from itself, or friction factor must be known al 
friction factor may be evaluated. 
It is evident that aging of the pipe tee. varying conditions of use site s0- 
salt a varying of absolute roughness which at present are not easily pre- 
Beep eed _ The experiments show that the surface of the pipe is changing with 
the age increasing the value of f. The factor e expressedininchesisa 


_ complete measure of the size, shape, and distribution of individual “et, 


’ The pipes ‘that are » rough in practice generally have no uniformity of the 
_ roughness element as assumed these by the authors of the equation. As there 
is at present no method for measuring or defining distribution of roughness - 
: element, the value of e is hypothetical.(1) In author’s examples illustrate — = 
steps i in an with the aid the modified di- 


x 


"instead = 5.64 


ge. 
0.456 


Sting 
Recommended use of Fig. 1 to find exponent in example 2m= s a 
°F agreement with f and e/D. It would be helpful to show symbols how to find 


_ The pipe resistance relations listed in the Table I as well as many other _ 
_ equations known by the writer, were based on experimental data with the flow 


of water, therefore, in common language, in development ipe 
must be used term water or liquid versus fluid. ( 
“A liquid is a body in that state in which the particles move re freely y among 
_ themselves, but remain in one mass, keeping the same volume but taking __ 
a always the form of the containing vessel; a liquid is an inelastic fluid; a gas 
is an elastic fluid that tends to expend to the utmost limits of the containing © Pi 
space. All liquids are fluids, but not all fluids are liquids; air and all the tl 


as a gases are fluids, but they are not liquids under ordinary circumstances, _ 
_ tho capable of being reduced to a liquid form by special means, as by cold > 


and Water at the ordinary temperature is once a figid and a 


AS 
, f, and e exponents from the modified erent in the ole 
eisbach are for designers | to | know these values 


The values of. e, 


8 
| wi 
a an 
an 
4 
ha 
4 te: 
= 
= 
s 
— 
= 
— 
— 
= 


min: it It: is ‘that, if accuracy alone were the de- 


ciding factor, any of the values given from the modified diagram, as wellas _ y 
City n, and exponents in other suggested formulas which are in —- use, — 


properly applied would be equally satisfactory. ba 


Victor Streeter. 


ASCE, Vol. 101 PP 682-704, 1936. 


q 948, Ed. 1942 


_ J. V. DAWSEY,! M. ASCE. —Professor has information 4 
_ which was heretofore scattered in many parts of the literature. Research to 


= methods for predicting or computing resistance encountered by = | 
- fluids flowing through a pipe is of interest in various fields of onguetting | 


Numerous to determine the relative roughness of a pipe, 


have been completed throughout the world in recent years. Some of these . 


tests were made ‘for the purpose of correlating the general resistance equa- 7 
tions with friction and roughness data-collected in model and prototype tests. 


=a on large conduits by forcing air at 1B velocities through . 
te-lined and unlined rock tunnels. Viparelli( 


The late Andre Jorissen(1) of Belgium experimented with relative roughness *< 


of Italy recently has 


The Waterways | Experiment Station is using a method devised by Mr. Frank 


‘duit linings. Two-in.-square plaster cast impressions are made at selected — a 


Campbell, (3) F. . ASCE, to determine the equivalent sand roughness of con- 


’ locations in the conduit that represent the roughness of some percentage of 
total conduit lining. Vertical-displacement measurements of the cast sur- 


An 1 Ames dial is presently used for surface-profile measurements on <i 
plaster casts after aninitial attempt to measure profiles electrically. The — 
point of contact on the shaft of an electrical contour analyzer was modified 


‘@ 4 = profile yield an average surface roughness of the conduit lining repre- _ - 


q ‘4 a, by using the point from a ball-point pen. This pen point indented the surface — : 


| ‘of the soft plaster casts under force of the analyzer probe. The Ames dial 


_ Shaft was also modified by the addition of a similar pen point, but did not in- 


dent the cast surface during measurements. For profile measurements, the 


} casts are firmly anchored in a mechanical vise and traversed under the shaft 


of the ; Ames dial. If the root- mean-square deviation of the surface “yy 2 


: - measurements is interpreted as the equivalent sand-grain roughness ¢€ (one 
a F half the diameter of the sand grain k in inches), then the friction factor f, at 
% the Reynolds number of the test, can be computed by using the general pipe im 


x Roughness casts of areas inside the 


= me Hydr. Engr. (Research), Prototype Section, Hydr. Analysis Branch, 7 
Hydraulics Div., U. Engr., Experiment Station, 


£ 
of the Enclish Lancuage _ 
experimented With a photographic determination Ol 
a 
a 
— 
— 


-square ‘impressions of the conduit with elastic dental - -impression 
material. Plaster casts were then reproduced from the impression material 

- negative. Equivalent roughness € for the conduit lining was determined from : 
root-mean-square deviations of the surface profile measurements on the 

rE = casts. _ This value of € was then used in the Colebrook-White equation, which © 
includes the effect of Reynolds number, to compute the Darcy-Weisbach fric- a4 


measurements pod roughness casts were recently at Enid 

| Dam. (5) The value of { determined from the casts was only 15 per cent a ae. 
than the value computed from the slope of the friction gradient. 


computed values of roughness, € €, the e corresponding friction factors differ _ 
- much less because of the exponential - logarithmic relationship between f and 
. | With s some additional refinment in sampling the conduit surface | and : 
analysing the plaster casts, it should be possible to ‘determine conduit rough- — 
accurately from plaster casts alone, including 
types of forms and construction techniques ar are used. 
Additional plaster - cast roughness data are needed to aid in the analysis 
= and interpretation of the equivalent sand-grain roughness relationship in | Sect 
: turbulent flow. Information from anyone making casts in conduits where the _ 
aie ‘friction gradient is known or can be measured would be appreciated. a 
Oe: ti It appears that parentheses were inadvertently omitted in the Colebrook- ~- 


a equation in Table I of Professor paper. - The equal equation should 


Jorissen, Andre, Contribution a L’ Etude des de Charge Continues 
dans les Conduites Circulaire, Bull. Contre Etude Constr. Genie Civil _ a 
2. Viparelli, M., Photographique Des Rugosites 
Determination of Roughness), Universita di Napoli, Recherches 


Hydrauli ues » 1957, - 250 , No. 15. 
Campbell, Frank B, Chief, Hydraulic Analysis Branch, 


“Pressure and Air Demand Tests in Flood- 


; Red River, Oklahoma and Texas.” MP No. 2-31, April 1953, B. Guyton, _ 
U.S. “Army Engineer Waterways Experiment Station, 
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4 — 


River, Mississippi. No. 2- 510, 1959, Cc. J. Huval, s. Arn 
s Exper 
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HYDRAULIC ANALYSIS OF SURGE BY COMPUTER 


important problems. 


The use of digital computers 


The use of digital computers for advanced research on surges in — 
_ systems of surge tanks is clearly demonstrated by the author. Digital com- ua 
_ puters have been used in hydraulic research by several authors especially - " vs 
| The paper by McCaig and Jonker, is a typical example of what can be | hay 

_ achieved with this new tool. The writer too has used a digital computer for — 
‘solving a simple but tedious problem. it was suggested that a small hydro — 

€ ‘power station could be provided with a surge tank too narrow to cope with 

large surges, provided the power sets were hand operated and put slowlyon 

R. the line, the operator observing very definite rules concerning the turbine a. ; 
+ gates opening times. _ There was no turbine governor and the load on the sta- 

was supposed to remain constant for long periods. 
ny When a large number of possible load combinations and possible errors i by ve 


operator were considered and effect on surges the 


: = de France. . The precision given by a direct calculation on the 

| digital computer goes far beyond of what can be achieved with a graphic — re 

_ diagram. In the case described by McCaig and Yonker such high precision —_ = 

_= required for comparison of cost estimates for several design alternatives. s. 
Paynter has been using the electric analog method for research on surges 
and on water hammer. His methods are highly suitable for very general re- 

F problems concerning concentrated pressure and head losses" (losses 

- occurring at this junction of a surge tank with the pressure tunnel) have to be 
analysed or when turbulence is involved, experiments on model tests is om 7 


‘The ‘Thoma. condition for su surge stability 


vo The e author correctly y remarks that ‘the whole si surge stability theory | estab- 


lished during aad Past 40 years by Thoma, Calame and Gaden, Escande and — 
a. Proc. Paper 1996, April, 1959, by Mr. Nicholas Barbarossa. phi a 


Co., Ltd. , Rugby, | England. 


7 


q 


| 
_ 
q a 
; 
7 
= 
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1959 
— is ee on the seoumstion that the demand for power varies abruptly, | 


} then stays absolutely constant for the period under observation and that he s 


4 
gg an basic assumption to be over- -simplifying. - According to Fig. 10, 11, and 12 ae 
_ _ the paper, the gates continue to open to “full gate opening” beyond the point P 
ae = should be remarked here that the points on Fig. 10 and 11 on which the — — 
discussion is based correspond to the very beginning of the surge: time _ 


A power generation could not be established during ‘the rest of the r period. De- 

pressed water levels will raise again. In addition, the governing equation re- 
@ lating speed variations to the moment of inertia of the > rotating masses could 
fi possibly be introduced asa an additional condition to. the computation: it m may 


a a show how the governor tends to stabilise the power output and additional re-_ 
search is suggested. main point to be stressed here is the authors be- 


lief shared by the writer, that further research on surge stability is definitely | 


5 - wanted. A first phase of research has been concluded: simplified mathemati- | 
ee. theories expanding ‘Thoma’s stability creterion to all types of surge > tanks 
7 and systems of surge tanks have been worked out. The Routh-Hurwitz mathe- — 
= matical stability criterion and the Liapunoff theorem have been used success- 
os fully by skilled mathematicians. All this approach work is based on the same 
— assumption of abrupt change in the load demand, which then remains | 


as: ‘How does this basic assumption fit in with the behaviour of real electric 


_ systems? — This is the question facing us for the second phase of the re cole 


search. A certain amount of preparatory work has already been done: 

_ Chevalier and Hug (Electricité de France) have tested the | Cordeac surge 

tank: They found that the governor was reacting correctly to the p power de- 


a - mand but with a time lag of about 10 : seconds. This time lag was found to a 


 Cuenod Gardel tested the differential surge tank 
Oelberg (Switzerland). , By modifying the actual surge tank they caused it to { 

7 become unstable. They tested different types of feed back linking surge J 
; and line voltage to the frequency in order to re-establish the stability of the | 


agi ‘ Finally Aillert (Electricité de France) has" designed: and | built several i | 


surge tanks, which are definitely unstable but where the surges are con- has J 
Ae These are definite efforts to approach the ‘Stability problem from a com- mp 
pletely new angle. More light ‘could be thrown on the way governors - react 


and surges can be controlled. Possible savings on the tank volume could be_ | 
achieved thanks to a closer of the physical in large elec-— 


It is suggested that a comprehensive paper summarising all these new 
2 4 “forts should analyse th the up to date position of the and 
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Ginn 


problems involving transients. A.S.M.E. ‘Paper No. | 58 

2. qneger, —— in surge tank design. Proc. Inst. Mech. Eng., a 


C. The double surge tank systems. Wi Power, Vol. 


August 1957, and La Houille Blanche No. 4, September 4 


Jaeger ontribution to the stability theory of systems of surge tanks. _ 
S.M.E. Paper No. 51 - A.65; Trans. A. A.S.M.E. Vol. 80, No. 7, — 


5. Chevalier, J. and Hug, M. ‘Tests on the Cordeac surge tank, concerning 

La Houille Blanche, Vol. 12, No. 6, reaorecees 1957. 


Thoma ‘Ss 


. Cuénod, M. and Gardel, A. Communication du Laboratoire @’ Hydraulique 4 > 


_ Ecole Polytechnique, Université de Lausanne, 1953, and i La Houille ~~ 
Blanche No. . 3, June 1954. 7 


‘ 
= Bouvard and Molbert. la Houille | , Blanche, Vol. 8, ‘No. 2, May 1 1953. | 


| 
10r has warmly Congratuia e or his successful use of a 
| ze computation and for his criticism of the assumptions on which __——— 
t 
— 
| 
| 
a 


oF that 


” into the 
yance factor Kg. There seems little to gain from including Manning’s | a ao 
coefficient in order to determine some nebulous quantity of conveyance which — : 
4 is of little use in itself. - discarded this concept some time ago(1, 2) 
: with the result that K(=1.486 AR 2/3) is determined solely by the geometry — 


of the channel c cross- }-section and the stage. Such values ‘when graphed o on na 


the discharge is given ai te: Q= = 


was included. — 

REFERENCES 


a +» Silvester “Flow in in Open Channels” J. _ Inst. - Engrs Aust. Vol. 25, +1953, 
* Silvester “New Graphs Aid | Calculations” ‘Water 
1954, $39. | 


= 


Proc. 1997, April, 1959, and K. B. ‘Schroeder. 
Senior Lecturer ‘in Civ. of 


— 
aie resuiled Straip eS ana uUrves depending upo id 
the shape of the channel. Graphs were supplied for trapezoidal, parabolic, § im im Sie 
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@ectangular, circular and U shaped channels. Once K has been so determined 
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